The effect of parallel static and microwave electric 
fields on excited hydrogen atoms 

D Richards 
Faculty of Mathematics 
Open University 
Walton Hall 
Milton Keynes, MKT 6AA. 

February 1, 2008 



Abstract 

Motivated by recent experiments we analyse the classical dynamics of a hy- 
drogen atom in parallel static and microwave electric fields. Using an appropriate 
representation and averaging approximations we show that resonant ionisation is 
controlled by a separatrix, and provide necessary conditions for a dynamical reso- 
nance to affect the ionisation probability. 

The position of the dynamical resonance is computed using a high-order per- 
turbation series, and estimate its radius of convergence. We show that the position 
of the dynamical resonance does not coincide precisely with the ionisation max- 
ima, and that the field switch-on time can dramatically affect the ionisation signal 
which, for long switch times, reflects the shape of an incipient homoclinic. Sim- 
ilarly, the resonance ionisation time can refiect the time-scale of the separatrix 
motion, which is therefore longer than conventional static field Stark ionisation. 
We explain why these effects should be observed in the quantum dynamics. 

PACs: 32.80.Rm, 33.40. -Hf, 34.10. -Hx, 05.45.Ac, 05.45.Mt 
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1 Introduction 



A strong electromagnetic fields can perturb an atom in many unexpected and compli- 
cated ways that are difficult to understand. If the atom is initially in an excited state 
usually a large number of unperturbed bound states are coupled, making the numer- 
ical solution of Schrodinger's equation difficult. Moreover, the corresponding classical 
dynamics is normally partially chaotic — meaning that there are both unstable and 
stable orbits close to the initial unperturbed torus — and the wave function mimics 
this behaviour, see for instance Leopold and Richards 1994 and Richards 1996b, thus 
making the interpretation of numerical solutions difficult. 

The investigation of the effects of strong periodic fields on an excited atom dates 
from the original experiments of Bayfield and Koch (1974) which showed that a rela- 
tively weak field could produce a multiphoton transition into the continuum, contrary to 
the received wisdom of quantal perturbation theory. The subsequent history of this in- 
teresting problem is told by Koch (1990). In 1974 conventional quantal theory required 
high-order perturbation theory to describe the 80 photon jumps to the continuum of 
this early experiment. This was, and remains, an impossible calculation, but theory was 
rescued by Delone et al (1978) who proposed that the classical ionisation mechanism 
involved diffusion of the electron through atomic states highly perturbed by the field. 
A year later Meerson et al (1979) proposed a different classical diffusion approximation. 
In the same year Leopold and Percival (1979) used a classical Monte-Carlo trajectory 
method to estimate classical ionisation probabilities and obtained qualitative agreement 
with experiment. In 1985 the experiment was repeated with better control of all impor- 
tant parameters and the comparison between these results and a classical Monte-Carlo 
simulation, van Leeuwen et al (1985), showed remarkable agreement. 

Since then our understanding of the dynamics of this type of system and the re- 
lationship between classical and quantal solutions in different parameter regimes has 
developed. For instance we now know that the scaled frequency — the ratio between 
the driving frequency and the Kepler frequency of the initial unperturbed motion — is 
one of the most significant parameters and that there are six separate scaled frequency 
regions in which the dynamics has quite different characteristics, see for instance Koch 
and van Leeuwen (1995). 

Linearly polarised fields with low scaled-frequencies were considered in Richards et al 
(1989) and there it was found that at particular scaled frequencies quantal effects were 
important due to resonances between two adiabatic states, see also Dando and Richards 
(1993). Low frecjuency elliptically polarised fields, Bellermann et al (1996, 1997) and 
Koch and Bellermann (2000), showed the existence of complicated resonance structures 
that, on the other hand, can be explained using classical dynamics, Richards (1997). 
When the scaled frequency is close to unity the main classical resonance island plays a 
dominant and similar role in both the classical and quantal dynamics, except at certain 
frequencies where scars produce differences, Leopold and Richards (1994) and Richards 
(1996b). At higher scaled frequencies classical dynamics fails, as predicted by Casati 
et al (1984), and demonstrated by Galvez et al (1988) for linearly polarised fields and, 
for elliptically polarised fields, by Wilson (2003). These different behaviours have been 
the subject of several reviews Jensen et al (1991), Koch (1990, 1995), Koch and van 
Leeuwen (1995). 

In this paper we consider the effect of strong, parallel static and oscillatory fields on 
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an excited hydrogen atom, the work being motivated by recent experiments of Professor 
Koch's group. These two fields affect the system in comphcated ways. Roughly, the 
static field splits the hydrogen degeneracy introducing another frequency to the sys- 
tem with which the external field can resonate. Classically this means that the Kepler 
ellipse moves periodically and this motion can resonate with the driving field. Such 
resonances can enhance the ionisation probability so a significant part of this paper is 
devoted to understanding the dynamics of these resonances and the mechanism caus- 
ing enhanced ionisation. Preliminary experimental results and some comparisons with 
classical Monte-Carlo calculations are given in Galvez et al (2000); a more detailed dis- 
cussion of the experimental method with more extensive results and comparisons with 
classical calculations is provided in Galvez et al (2004), subsequently referred to as I. 

In these experiments excited hydrogen atoms are subjected to strong parallel static 
and periodic fields, so the external force on the electron is {Fg+Ff^ cos J7t)z: the accurate 
theoretical description of this system present a challenge. The main observed effect of 
this field combination is to produce an ionisation probability that, for fixed J)i and f2, 
rises steadily as Fg is increased, but which is punctuated by a series of approximately 
equally spaced sharp local maxima, see for instance Galvez et al (2000, figure 1) and 
figure n of this paper. Each local maximum is produced by a resonance between the 
periodic part of the perturbation, zFf^ cos ilt, and the mean rotational motion of the Ke- 
pler ellipse induced by the static part, zFs. These peaks are relatively narrow but their 
shapes, widths and heights vary significantly with the system parameters. Moreover, 
they disappear at certain values of Ff^ and Fg, for reason that are sketched in Galvez 
et al (2000); similar explanations are provided by Oks and Uzer (2000) and Ostrovsky 
and Horsdal-Pedersen (2003). 

Since publication of the original Galvez et al (2000) paper we have made detailed 
comparisons between the ionisation probabilities of the experiment and classical Monte- 
Carlo simulations, which will be published in Galvez et al (2004). These comparisons 
show remarkable agreement in many instances, but also some differences, some of which 
may be attributed to unquantifiable differences between the assumed substate distribu- 
tion. In experiment and calculations we observe systematic differences between both the 
positions of the resonances and of their disappearances and the predictions of the simple 
theories. In an attempt to understand these differences more detailed numerical inves- 
tigations of the classical dynamics were undertaken and these calculations show that 
the dynamics underlying the apparently simple, averaged results seen in figure^ is, in 
fact, very complicated: even an accurate theoretical prediction of positions of the local 
maxima in the ionisation probability is fraught with difficulties in classical and quan- 
tum dynamics. Further, it transpires that the nature of the resonance that causes these 
maxima is unusual in that the classical resonant island moves as the field is switched 
on: in some circumstances this means that the classical ionisation probabilities reflect 
the development of a homoclinic tangle, see figure [TBI and [T^ The size of the quantum 
numbers needed for the quantum dynamics to mimic this behaviour is not known. 

This paper has two main aims. First, to provide a theoretical understanding of 
the classical dynamics and to isolate those features that determine the nature of these 
resonances, for example their positions and temporal development. Second we need to 
derive an approximate Hamiltonian that leads to a numerically tractable Schrodinger 
equation. Both aims are achieved by using a representation in which coupling between 
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basis states is relatively small. 

In section |5] notation is defined and we present some numerical results illustrating 
the theoretical problems that need to be solved. In section 13 the approximations are 
developed; this is algebraically complicated but, in essence, it is simply a standard use 
of perturbation theory and averaging approximations. However, for two reasons care is 
need with this theory. First, the convergence of the perturbation expansion needs to be 
understood because the fields used are strong and, it transpires, sometimes beyond the 
radius of convergence of the series. Second, it is shown that particular observable effects 
are produced by subsets of terms and it is necessary to determine their origin because 
only some subsets can be computed to high order — this means that some effects, 
such as the resonance positions, can be computed to high order, but others cannot. 
This analysis is necessary partly to extend the earlier approaches beyond first order 
and partly to show how ionisation occurs, a feature missing from all earlier theories of 
this system. Furthermore, from this analysis emerges an approximate Hamiltonian that 
should facilitate a quantal calculation. 

In section 0] the mechanism connecting the dynamical resonances, described in sec- 
tion 13 with ionisation is described. A significant result of this analysis is that the 
positions of the dynamical resonance is not precisely the same as the positions of the 
local maxima seen in ionisation curves; although small, an estimate of the difference 
seems beyond current theory. In this section we determine some necessary conditions 
for a dynamical resonance to affect ionisation and also obtain accurate estimates for 
the position at which the resonances observed by Galvez et al (2000) disappear, see also 
Schultz (2003) and Schultz et al (2004). We also discuss resonance widths and show that 
for the current experiments, which involve an average over substates, the width is not 
due to the variations of the resonance position with substate (which is relatively small), 
but is due mainly to non-adiabatic dynamical effects, that are difficult to quantify. 

In section 121 we analyse the time required for a resonance to develop and show this 
to be relatively long. Further, we observed that the nature of the resonance island 
means that the classical ionisation probability can be significantly affected by the field 
envelope. In particular in some circumstances the ionisation probability can reflect an 
incipient homoclinic tangle that develops as the initial state moves slowly through a 
separatrix. 

2 Notation 

The Hamiltonian for a hydrogen atom in parallel static and microwave electric fields, as 
in the experiments of Galvez et al (2000, 2004) (the latter henceforth being referred to 
as I), is derived in Leopold and Richards (1991) and, provided the field envelope \{t) 
changes sufficiently slowly, is given by 

1 

H = —Yi' ~ - + F{t)z, F{t) ^{F, + Ff, cos nt)X{t), (1) 

where /i is the atomic reduced mass, e the electron charge and X{t) is the envelope 
function describing the passage of the atom through the cavity. This Hamiltonian has 
azimuthal symmetry so the z-component of angular momentum, /„i, is conserved. For 
the particular experiments described in I, X(t) has the 16-113-16 configuration, meaning 
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that it rises monotonically from zero to unity in 16 field periods, remains constant for 
113 periods and then decreases monotonicaUy to zero in 16 periods. In aU calculations 
reported here the initial rise over Na field periods is taken to be 

X{t) = x^{2- x^), t^^^^a;, < a: < 1, 

and the decrease as the field is switched off has the same shape. Classical ionisation 
probabilities are normally insensitive to small changes in X{t)] exceptions to this rule 
are discussed in sectional 

For excited atoms it is convenient to use units defined by the initial unperturbed 
Kepler ellipse which has semi-major axis a = /Q/(/ie^) and frequency uik — Me'*//g, 
^0 = noh, where Iq and no are the initial values of principal action and quantum number 
respectively: scaled units are convenient because the magnitude of most scaled param- 
eters that produce similar physical effects change little with uq. The scaled frequency 
and field amplitude are defined by 

f7o = — = ^ = (0.00533757no)' , 

ujK Aie-* GHz 

Fo = ^ = ^ = -— (0.00373535no)* . 
e /i^e° V/cm 

The scaled time is to ~ ujut and a scaled action / is I/Iq. In the current experiment 
n = 8.105 GHz, so rio = (0.010722no)^. In the following we use the symbols Fg and 
for both the scaled and actual field amplitudes and t for scaled and actual time: this 
misuse of notation avoids a clutter of subscripts but should not cause confusion because 
scaled quantities are dimensionless. 

2.1 Some numerical results 

Before describing the theory we show the results of a few classical calculations in order to 
provide the reader with an idea of the features that a theory needs to describe. For the 
present calculations a Monte-Carlo method, as described in Abrines and Percival (1966), 
is used in which N initial conditions are chosen from a microcanonical ensemble: if M 
of these orbits ionise the estimate of the ionisation probability is Pi = M /N. Without 
stratification the standard deviation of this estimate is (Hammersley and Handscomb, 
1964) 

N ' ' N' 

meaning that there is a 68% and 95% probability of the true result being in the ranges 
{Pi— a, Pi+a) and {Pi — 2a, Pi+2a) respectively. In the present calculations the sample of 
initial conditions is stratified by dividing the range of each variable into equally probable 
intervals and choosing, from a microcanonical distribution, one point at random in 
each sub-interval, see Abrines and Percival (1966). Stratification reduces the statistical 
errors and sample calculations suggest that with this form of stratification the true 
value of a is about half the above estimate. For the numerical integration of Hamilton's 
equations the problem associated with the Coulomb singularity is circumvented by using 
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the regularisation method described in (Rath and Richards 1988): numerical integration 
was performed with the NAG routine D02CAF. 

In the first figure is shown a 'typical' classical ionisation curve in which F^^ and fio 
are fixed and the variation of the ionisation probability, Pi, with the static field Fg is 
shown: here a microcanonical distribution of substates is used. For this illustration we 
choose rJo = 0.0980, (no = 43) and F^ = 0.10: the field envelope was 16-80-16 and 1296 
orbits, for each value of Fg, were used. 




0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 



Figure 1 Ionisation curve for SIq = 0.0980, (no = 43), = 0.1. 

The broad features are clear: Pi = for Fg < 0.013 and Pi = 1 for large Fg, typically 
Fg > 0.2. As Fg increases between these values the steady increase in Pi{Fg) is punctu- 
ated only by sharp local maxima at almost equal intervals in Fg, at Fg = 0.0316, 0.0620, 
0.0916 and 0.1203; a close inspection of the data shows another small maximum at 
Fg ~ 0.148, marked by the arrow. The small amplitude undulations seen in this ion- 
isation curve are assumed to be caused by the statistical errors mentioned above, and 
provide a visual estimate of the magnitude of these errors. The variation in the scale- 
length of these oscillations in Fg is due to variations in i<"s-interval used in the calculations 
which was smallest near the local maxima of Pi. 

We denote the positions of these maxima in Pi{Fg) by J-'i'''' . Theory associates these 
maxima with a resonance between the driving field and the precession of the atomic 
Kepler-ellipse and we denote the position of these resonances by Fg . The two field 
values Fg''^ and jFi"*^ are approximately the same but are not identical, which is why it 
is helpful to use different symbols. 

Elementary consideration, section 1X51 show that Fg''^ ~ (in scaled units), 

giving 0.0327, 0.0653, 0.0980 and 0.131, so the relative differences between rioj/S and 
are 3.5%, 5.3%, 7% and 9% respectively. The more accurate theory developed in 
section lO gives F^'^ = 0.0316, 0.0623, 0.0904 and 0.112, for j = 1-4, with relative 
differences of 0.5%, 1.3% and 7%, respectively, for j = 2, 3 and 4, with the values for 
J = 1 agreeing to three figure accuracy. 

The averaged ionisation probabilities disguises a richer and more complicated be- 
haviour. In the next three figures are shown ionisation probabilities from a given sub- 
state, I2 — Im = 0.2 — these scaled actions are defined below, equation |S1 and |^ Here 
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flo = 0.0528, {no — 35), and 1600 orbits, for each F^, were used. The arrows point to 
where F^+^ 



Fs = f!:^} — F,, where F^:^} is defined in equation 1211 befow. 



F =0.14 




F =0.15 
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Figure 2 lonisation probabilities for Qq = 0.0528, ^2(0) 
16-50-16 envelope and for various values of F/j,. 



0.2, (/e(0) = -0.4), with the 



At the lowest microwave field, F^^ = 0.13, a number of local maxima are seen: those 
labelled j = 1-4 can be associated with a dynamical resonance, and most are visible 
for Ff^ — 0.14 and 0.15. For = 0.13 there are also a number of other local maxima, 
at Fs — 0.0430, 0.0535 and 0.0705 the origin of which is not known, but in section [S] 
results are presented which suggest that the maximum at Fs ~ 0.0430 is a non-integer 
resonance, equivalent to j = 2|; In this example Pi = for Fg < 0.016 and Pi = 1 for 
Fs > 0.076, with an underlying steady increase in Pi{Fs) for Fg > 0.045. 

For = 0.14 the four labelled maxima of the previous figure persist but have shifted 
slightly and the j — 3 maxima has split: we have confirmed that the latter effect is not 
due to statistical sampling errors. In other calculations with the shorter rise time of four 
field periods there is no split, suggesting that it is caused by the field envelope: this and 
other effects of the field envelope are discussed in the sectional The j = 1 maximum 
near F^ = 0.0166 is lower and the theory developed below shows that it disappears 
completely for F^ ~ 0.147, Fg ~ 0.0161. A new feature in this graph is the large value 
of Pi at Fg — 0, with Pi falling rapidly to zero as Fg increases to 0.007. For Fg > 0.035 
the underlying trend in Pi is a steady increase to unity, but all other structure is not 
understood. 

For = 0.150 the three labelled maxima persist, but have shifted and broad- 
ened. Now Pi is large for Fg < 0.013 and Fg > 0.04 and there is a new maximum at 
Fs = 0.0192, the existence and magnitude of which depends upon the switch-on time, see 
figure 01 we have confirmed that, at these low frequencies. Pi is affected insignificantly 
by the switch-ofF time, provided the total interaction time is sufficiently long. 

The numerical details of the resonance positions seen in these figures are given in 
tabled Here the values of J^i"'-' are computed using a grid AFs = 0.0001 and defining 
^i"''' to be the value of Fg at which Pi{Fg) has a local maximum. If Pi > 0.99 for a 
range of fields, the range is quoted and sometimes there is more than one clear maximum. 
For fixed F^ and increasing j there are clear systematic differences between these two 
values. It is not known what causes these differences; the breakdown of the perturbation 
expansion used to compute Fg'''^ maybe significant, but also the discussions in sections f4.2l 
and El show that the relation between the values of Fs"*^ and ^i"*^ is far from simple. 
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Table 1: Values of J-s , computed by the method described in the text, and F, given by 
the theory described in section 
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F^ = 0.15 
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Monte Carlo 
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Monte Carlo 


theory 


Monte Carlo 
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1 


0.0166 


0.0166 


0.0162 


0.0164 


0.0158 


0.0160 


2 


0.0327 


0.0330 


0.0320-0.0323 


0.0323 


0.0314-0.0316 


0.0314 


3 


0.0486 


0.0487 


0.0464,0.0483 


0.0474 


0.0455-0.0483 


0.0455 


4 


0.0645 


0.0631 


0.0632-0.0641 


0.0607 







The regular but not equally spaced series of four minima beyond the arrows, in figure |21 
have, at present, no dynamical explanation. Numerical evidence, however suggests that 
some of this structure is caused by the field switch, with shorter switch times removing 
most of the structure and longer times producing more, see figureOl the duration of the 
centre part of the envelope is irrelevant provided it is long enough, see sectional 

In the following two figures we show ionisation curves for the same parameters used 
in the right hand panel of figure [21 but with difi^erent envelopes. 
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Figure 3 Ionisation probabilities for Ho = 0.0528, /2(0) = = 0.2, (/^(O) 
Fjj, = 0.15 and various envelopes. 



-0.4), with 



For the shorter switch, 4-50-4, Pi 1 for Fs < 0.014 and F^ > 0.052 and there are local 
maxima only at the j — 1 and 2 resonances where J^s^^ = 0.0158 and J-i^'' — 0.0315. 
For the longer switch, 40-50-40, the above two maxima persist, at the same field values, 
but now the ionisation curve shows a great deal of other structure, some of which is 
sensitive to the switch time. We discuss other effects of the field switch in section |H1 



3 Theory 

In the situations of interest here the field frequency is small by comparison with the 
Kepler frequency, so the scaled frequency Op, equation [3 is small; typically it varies 
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between 0.04 and 0.11 for no between 32 and 45 for the 8.105 GHz cavity. Hence the 

variation of the field, F(t), is slow by comparison with the electron orbital motion and 
an averaging approximation may be used to remove one degree of freedom, but first it 
is necessary to choose the correct representation. 

The fields encountered here are sufficiently large to couple together many states of 
the field free atom and to ionise some of these states, so useful theoretical descriptions of 
the experimentally observed signal must include a mechanism for ionisation and, ideally, 
should use a representation in which coupling between bound states is relatively small. 

In order to understand the magnitude of this problem, and to motivate the following 
analysis, we consider the static field ionisation of the one-dimensional atom with the 
Hamiltonian ^ 

Hi = _ ^ + x{t)Fz, z>0, 
2fj, z 

where F is constant and here A increases monotonically from zero to unity over a time 
long compared to a Kepler period. We denote the quasi-bound states of Hi when A = 1 
by |nF). If initially the atom is in the state |no 0) there is no classical ionisation pro- 
vided F < Fcrit, where, in scaled units, F^nt = 2'^"/{3tt)'^ ~ 0.1298, Richards (1987), 
and complete ionisation if > i^crit- In quantum mechanics tunnelling decreases these 
thresholds by an amount that depends upon the interaction time and the initial prin- 
cipal quantum number. For an excited Id hydrogen atom initially in state \nF) the 
probability of remaining bound at time t can be deduced from the relations derived in 
Richards (1987) and behaves approximately as, 

Pfc(n,t) ~ exp(-R), n^r ^ i- exp (^~2.58nc(n) . < ^crit) 

where for n > 20, c{n) ~ 1/(1 + 1.65/n + 173. l/n^- 249. 5/r?;'') is derived bv a numerical 
fit to the theory. This probability approaches a step function as n — > oo. For t = 2?™^, 
that is one Kepler period, the probability Pb decreases from 0.9 to 0.1 as F increases 
over an interval 7Fcrit where 7 = 0.2, 0.1, 0.04 and 0.02 for n = 5, 10, 30 and 50, 
respectively. 

For F < Fcrit the state jng F) can be approximated by a linear combination of the 
field free states; the matrix elements (no F \ nO) provide some idea of how many unper- 
turbed states are required to accurately describe the wave function in the presence of a 
strong field. This matrix element is estimated by Richards et al (1989) where it is shown 
to be significant for no < n < m, where m ~ no 

(2F)-i/4 (= i.5„o for F = 0.1). Thus 
any realistic approximation using a basis of unperturbed states requires about 2no states 
in a Id system and about ng states for each of the n values of the azimuthal quantum 
number, m. The method used by Robicheaxix et al (2002) avoids these problems, but as 
no increases the number of grid points increases and the computational time increases 
commensurately. 

The theory presented here minimises coupling between basis states by describing the 
motion in a basis that diagonalises the static-field, or Stark, Hamiltonian 

1 

Hs = — \-Fz, F = constant (> 0). (3) 

If |nF) is a bound eigenstate of Hs with energy E^{F) then a basis that may be used 
to describe the bound motion of the time-dependent Hamiltonian is obtained simply by 
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replacing the constant F by F{t)] in this basis coupling between states is caused only by 
the rate of change F{t) not by its magnitude. For the examples of interest here it will 
be seen that this coupling is, in scaled units, 0{floF^), which is typically an order of 
magnitude smaller than the coupling between the unperturbed states. This method was 
first used in the present context by Richards (1987) and Richards et al (1989), where 
it was applied to the one-dimensional hydrogen atom and shown to explain important 
features of the three-dimensional experimental results. 

Because this approximation uses a bound-state basis the continuum needs to be in- 
troduced as an extra approximation, described later. In addition, the eigenfunctions 
(r I XI F) are not conveniently represented by simple functions, consequently approxi- 
mations to these are necessary. Finally, since canonical transformations are easier to 
handle than their corresponding quantal unitary transformations, it is easier to develop 
this approximation using classical dynamics and to quantise the resulting Hamiltonian, 
rather than tackle the quantum mechanics directly. 



The classical Stark effect 

The first goal is to find a suitable approximation to the generating function, S'(I, r, _F), 
for the canonical transformation to the angle-action variables of Hs] we also need ex- 
pressions for Hs and dS/dF in terms of these variables. This is a relatively routine, 
but complicated, calculation because it is necessary to expand to high orders in F . The 
main result of these calculations is the adiabatic Hamiltonian, defined by equation 1181 
below, which forms the basis of further approximations. 

The theory starts with the Coulomb-Stark Hamiltonian, equation |21 in which the 
force on the electron is static and in the negative z-direction. This Hamiltonian is 
separable in the parabolic coordinates; following Born (1960, section 35) we use the 
coordinates 

x = Ct?cos(/), y = ^'qsuY(t), z = ]^{^'^ -rf), £,>0,r]>0, 
sometimes named squared parabolic coordinates, giving 

In the following we assume F > and E < 0. The Hamilton- Jacobi equation is 



^iF{e-vn~^t^e' = 2fiE{e+v'), (5) 



the general solution of which defines the generating function ^(I, r, F) for the canonical 
transformation to the required angle-action variables. This equation is separable so 
S = Si (C) -I- 5*2 (?7) + /m0 where 

SiiO = lf{-f^Fe-2f,m' + 2a,f,e'e^l'y, (6) 
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where ai and ai are the dimensionless separation constants that satisfy a\-\- ai = 2, 
with ai > and 0L2 > 0. Motion in the ^-direction can be bound or unbound whereas 
motion in the 77-direction is always bound. The bound motion is restricted to the regions 
< ^1 < ^ < ^2 and < r/i < ?7 < 7/2 where ^1 and jji are zero only if Im = 0. 
The action variables are defined by the integrals 

h = - f {-nFe -2^\E\e + 2a,i,^e - II) , (8) 

I, = ^ ^^{^,Fn'-2^,\E\r,' + 2a,^eW-Il)'^\ (9) 

and satisfy the relation Ii + 12 + \Im\ = In and Q < < In — \Im\- They are related 
to the usual quantum numbers n\ and n2 by 

Ik = {^k + 1/2)% with m + n2 + |m| + 1 = n, < n/j < n — |m| — 1. 

These equations relate {Ii,l2) to {E,ai) and may be inverted to give E and ai in 
terms of (/i,/2). However, for F ^ the integrals cannot be evaluated in closed form. 
One method of inverting these equations is to invert the series obtained by expanding 
as a power series in F. A method of performing these calculations is outlined in the 
appendix; the resulting algebra is complicated and performed using Maple. For reasons 
that will soon become apparent, we have computed these series to 0{F^'^), but here 
quote lower order expansions. The resulting perturbation series for the energy is 

4 00 

^(I) = -^+E^'^W^'' iln = ll+h + \lm\), (10) 
fe=l 



where 






Ei{I) = 












Esil) = 




Eiil) = 


rlO p4 

-102Ve20 (1829/^ - 1134/^/^ - 183/^ + (602/^ - 


-3787^)7,2+49/4) 


Esil) = 


"Hps + + ^^^^^ + + 


2207^)7f-2174), 


Ee{I) = 


-819^^11^32 {54726276 - 429903/^7^ - 16200/X 


- 69517^ 



+ (68515274 - 254707^7^ - 364507*) + (3907^ + 7657^) 7* - 3727^} , 
E7{1) = -^^^^^-^{7071885/6 - 15305617*7^ + 94915/2/;^ + 559377^ 



+ (15022837* + 214107^7^ + 661157* ) 7^ + (19477^ - 63217^) 7* + 9577^} , 

and 7e = I2 ~ Ii] this last action variable is related to the electric quantum number, 
7e = rieli, though in some text, for instance Bethe and Salpeter (1957), the electric 
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quantum number is defined to be ni — rt2, because coordinates are chosen so the force 
on the electron is in positive ^-direction; the value of le is the projection of the Runge- 
Lenze vector along Oz. Up to 0{F^) the above expansion agrees with the series given 
in Damburg and Kolosov (1983, page 45), see also Silverstone (1978), as ?i — + 0. Note 
that the odd components i?2fe+i(I) have a term linear in Ie\ it will be shown that these 
components determine the resonance position. The separation constant is, to 0{F^), 



1 

128M^ei8 



{ (171/2 - 15/f ) {II - II) - II (82/2 ^50^2 ^ 25/2„) } + OiF^) 



with a2 — 2 — ai. 

The angle- variables corresponding to the action variables defined in equations|Hland|ni 
are defined by 9k = dS/dlk, fc = 1, 2. It is shown in the appendix that the two angle 
variables can be expressed in terms of relations 

01 = V^ + Pi(^) + Qi(x), ^?2 = x + ^2(V^) + Q2(x), (12) 

where the two auxiliary angles, {tp, x) ^-re defined by the relations 

e - 1 +e,)-l {e, - c?) cos^, ^ i (^2 ^ ^^2^ 1 (^2 _ ^2^ ^^^^^ (^3) 

and where {Pkix),Qkix)) are odd, 27r-periodic functions with zero mean value; expres- 
sions for these, accurate up to 0{F), are given in the appendix, equation 1501 However, 
for reasons discussed below, we require these functions only in the limit F = and then 
we have 

6*1 = - CTi sin?/; - (72 sinx, ^2 = X - Ci sin - (72 sinx, (F = 0), (14) 



where ak = \//fc(/fc + /m)//n, A: = 1, 2. 
Dynamic Stark effect 

When the field amplitude F varies with time the function, S{I,r, F{t)), generates a 
time-dependent canonical transformation and the Hamiltonian becomes 

K.Eil.FmH§% (15) 

The first term of this is the Stark Hamiltonian, equation IIUI the second term is more 
difficult to find, but it is important because only this term mixes states. It is shown 
in the appendix that the function dS/dF can be expressed as a Fourier series of the 
following form 

dS °° 

— = ^ sin fcf/; + -Bfe sin fcx (16) 

k=l 

where the angles (V'lX) s-^e defined in equation 1131 and where the coefficients {Ak,Bfi) 
are functions of the action variables and F. It is important to note that there is no 
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term independent of both ip and x- In our applications Fit) = \{t){Fs + Ff^cosftt) 
and A/A il, so F ~ — AilF^sinfit and since, in scaled units il, Fs and are 
numerically similar and small, a second order approximation is obtained by evaluating 
the derivative dS/dF at = 0, which considerably simplifies the analysis: the following 
result is derived in the appendix, 

where 

G = (3/2 +I1+ 2/,„)cri sin ?/; - (3/i +I2 + 2I„,)a2 sin x - y {aj sin 2?/' - sin 2x) . 
It is more convenient to use new angle-action variables. 

In = h + I2+ Im , Ol = <t>n - 0e, 

le^h-h, 62^(Pn+(Pe, (17) 

im — Jrrii — ^ ^ra , 

SO, when F = equation 1141 gives 2(j)e — O2 — 9i = x~ i^- This gives an approximation 
we name the adiabatic Hamiltonian, 

dF 

K = E{I,F) + ^-^—G{InJi,l2,i>,x), (18) 

where E{1,F) is the Stark energy given in equation 1101 The angles ip x not 
conjugate to the action variables and to develop further approximations it is necessary 
to express all quantities involving ip and x terms of 61 and 62, using equations 1141 
This is most easily achieved by expressing each function as a multiple Fourier series, 

^ ' Si= — OOS2 = — 00 \ ^Sl S2 / 

where cifla = si^^si and 
fc 

*7r'^s2(so'2) [Jsi+k{scri) + Jsi_fc(scri)] , S = Si + S2 7^ 0, 
2s 

±^^, s = 0, si = ±1, and fc = 1, 
0, otherwise. 

The adiabatic Hamiltonian 1 1 81 is useful because the coupling term is O(ilo^)i) rather 
than 0(F) as in the original Hamiltonian, so the resulting Schrodinger equation may 
be solved using a far smaller basis. 

Hamilton's equations in the original representation are singular at r = and this is 
dealt with using a regularisation method, see Szebehely (1967) for a general introduction 
and Rath and Richards (1988) for an application to the perturbed hydrogen atom. The 
equivalent singularity in the adiabatic Hamiltonian occurs when = and (Ti + (T2 = 1 
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and this also needs to be removed. For instance the equation for ip is ip ~ [ip, K] and 
the right hand side is proportional to 1/J where J = (1 — cri cos'0 — o'2 cosx), which 
can be zero. A method of avoiding numerical problems when J is small is to define a 
new time, r, by the equation dt/dr = J, to give the equations of motion 



'^^^ /„K^ (^(1 - (72 cos x)G^ + ai cos ipGx) ^ 

/„K^- ^(T2 COS X + (1 — (71 COS V") COS V' G^ , 

di! , dK fdK dK\ ^ 



dr 

dh _ . dF 
d7 



dx r dK fdK dK\ ^ 

Tr = '"a77+la77"a77j'""^'°'^ 

dF /2/i+/,„ 2/2 + /™ . \ „ 

where k = I^/{2^^e^) and = /i + /2 + ^m- In figure |S1 below we compare ionisation 
probabilities computed using these equations and the original Hamiltonian, but first it 
is necessary to re-introduce ionisation into this approximation. 



3.1 Ionisation 

The adiabatic Hamiltonian, equation 1181 does not allow for ionisation because angle- 
action variables exist only for bound orbits. Ionisation therefore has to be included as 
an extra approximation which is described here. 

For static fields each classical state, or torus, labelled by the actions (/„, Ie,Im), has 
a critical field Fcrit such that it exists only if < F < i^crit: the approximation to i^crit 
given by Banks and Leopold (1978) is used here. Note that if = bound orbits 
exist for all F, Howard (1995), but those orbits that exist for large F are so special 
that they do not affect the current problem. Adiabatic invariance suggests that this 
behaviour persists for sufficiently slowly varying fields, that is small scaled frequencies. 
The extreme values of Fcrit occur when the atom is aligned along the field, /,„ — 0: 
in scaled units min(Fcrit) — 0.13 for (/e,/m) = {ln,0), and max(Fcrit) — 0.38 for 
(Je, Im) = i^In, 0). The variation of i^crit with le = I2 — Ii, for various values of /„, is 
shown in the following figure. 
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Figure 4 The scaled classical critical field P'crit as a function of 
the scaled action le = I2 — Ii, for various Im. 



In the adiabatic limit /„ and /g are almost constant, so we may define the time- 
dependent critical field i^crit(0 — Fciit{Ie{t) ^ Im) and assume that ionisation occurs 
at the time when the actual field F{t) exceeds this. If F{t) changes sign, as happens if 
> the quantisation axis changes direction and the relevant critical field is given 
by i^crit(— -^e, —Im) {F < 0) SO the combined ionisation criterion is: 



This approximation is accurate when the field changes very little during one Kepler 
period, that is JIq ^ 1, and then this criteria may be used to include ionisation in the 
adiabatic equations of motion, equations 1201 

A useful guide to the behaviour of the system is obtained by putting /„ (t) and /g (t) 
equal to their initial values. This gives two boundaries beyond which Pi — 1: 



I: Fs > F^ n{Fs) = 1 if > F+t(/e, /„) - F^ 

II: Fs < F^ n{Fs) = 1 if Fs<Ff,- F-,,{Ie,Im) or F^ > F+,^{h,Im) - F^ 



where /e = -^e(O) is the initial value of this action variable. 

For the ionisation curves shown in figure |S1 /e(0) = —0.6 and = giving F'^^-^^ — 
0.219 and F^^j^ = 0.142. Thus, since F^ ^ 0.15 the condition II gives Pi{Fs) = 1 if 
Fs < 0.008 or Fs > 0.069: these boundaries are shown by the arrows in the figures and 
are consistent with the calculations. Similar boundary-arrows are included in figures |21 
andini 

Non-adiabatic dynamics affect this simple picture in two ways. First, they blur and 
slightly shift the boundaries. Second, and more important, isolated resonances produce 
large changes in Ie{t) and can enhance ionisation at particular combinations of {Fs, F^), 
other than those defined by I and II above. These dynamical effects produce the peaks 
labelled j = 1-4 seen in figure El 

In the following figure we compare values of Pi computed using the original Hamil- 
tonian ^ the solid lines, and the adiabatic equations EUJ the dashed lines. Here 
fig = 0.0528 (no = 35), I2 = 0.2, = and J), = 0.15 and the arrows denote 



F{t) > F^+> =F„,t{Ie{t),I„^), ff 
F{t) < -F(T^)=-F„it(-/e(i),-/™), ff 



F{t) > 0, 
F{t) < 0. 



(21) 
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the values of — F^^^^ and F^^^ ^ ■ ^^'^ ^^^^ hand panel the field envelope is 4-50-4 
and in the right panel it is 16-50-16. 




0.02 0.04 0.06 0.08 0.02 0.04 0.06 0.08 



Figure 5 Comparison of ionisation probabilities computed using exact dynamics 
(solid line) and the adiabatic equation I2UI (dashed line) for two envelopes. See 
the text for the explanation of the arrows, which point to the borders outside of 
which Pi = 1, in the adiabatic limit. 

These figures show broad agreement between the two calculations, but there are three 
marked differences; consider the left hand panel. 

• The maxima j = 2, 3 and 4 are at different values of F^. This is because only 
terms upto 0{F'^) were included in the expansion of E{T), eauation llOl and is not 
an inherent inaccuracy of the adiabatic approximation. 

• The maximum j = 1 is not present in the adiabatic calculation because with the 
approximations used it disappears when Fg — flo/3 — 0.176 and i^^ = floj[j./3, 
k = 1, 2, where jj^, are the positive, real roots of Jj{x) = 0; with fc = 3 
this gives i^^t = 0.150. For the full Hamiltonian this resonance disappears at 
Fs ~ 0.0161, F^, ~ 0.147. 

• The main difference is the shift in the shoulder near Fg = 0.01. 

Apart from these differences the agreement between the two approximations is good. 
The same remarks apply to the right hand panel but now we see that the new maxima 
at Fs ~ 0.039 and some of the structure at Fg ~ 0.07 are reproduced in the adiabatic 
approximation. The adiabatic approximation also has a local minimum at Fg = 0.0025, 
not present in the 'exact' probabilities: however, similar behaviour is seen in the exact 
results for other parameters values, see for instance the right hand panel of figure |31 

These results, and other comparisons that cannot be shown here, suggest that the 
adiabatic Hamiltonian provides a good approximation to the true dynamics. This is 
important because, for the principal quantum numbers used in current experiments, the 
numerical solution of the Schrodinger equation derived from this Hamiltonian is a feasi- 
ble computational task, unlike that derived from the exact Hamiltonian using either an 
unperturbed or a static-Stark basis. In the quantal approximation ionisation is included 
by adding a complex part to the energies, which may be computed semiclassically, as in 
Leopold and Richards (1991) and Sauer et al (1992). 
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3.2 Averaged equations of motion 

The adiabatic Hamiltonian 1 1 81 can be simplified further by noting that the two natural 
frequencies of the motion are quite different, 

so I We I <C ti-'n — ll!k\ this means that the orbital elements of the Kepler-ellipse change 
relatively little during one Kepler period. Hence the first averaged approximation is 
obtained by averaging over </>„, which is most easily achieved by ignoring all terms 
containing (pn — {6i + 62) /"^ in the Fourier series ^1 This gives 

(sin-;/') — — — (72 sin2(/)e, (sinx) = — (Tisin2(/)e and (sinfc-!/;) — (sinfcx) =0, k>2. 

Substituting these mean values into the adiabatic Hamiltonian ll 81 gives the mean motion 
Hamiltonian 

Km = E{h,F{t)) ~ --^_A(/e)i?(/e)sin20e (22) 

where 

A(I,f - (/„ + l^f - II B{hf - (/„ - I„,f - II 

In quantum mechanics this approximation corresponds to ignoring all transitions be- 
tween states of different n-manifolds of the adiabatic basis. 

Numerical integration of the equations of motion derived from 1221 is not straight- 
forward because of the square root singularity of B{Ie) at le = ±(/„ — /,„). It is best 
accomplished by introducing the three-dimensional vector 

Z= (B(/e)cOs20e,S(/e)sin20e,/e), 

the components of which satisfy the commutation relations [Zi,Zj] ~ 2eijkZk and 
|Zp ^constant so the vector Z moves on the surface of a sphere. The equations of 
motion, Zk = [Zk,H], become 

where k = /,^/(2/i^e^). These equations are trivially solved numerically. 



3.3 The Resonance Hamiltonian 

The mean motion Hamiltonian [221 needs to be rearranged in order to extract a clear 
picture of the dynamics. Observe that the odd terms in the series for the Stark energy, 
equation 1101 contain components that are linear in /g and that these terms produce 
a slow secular change in 0e which physically corresponds to a rotation of the angular 
momentum vector about the Runge-lenze vector. We shall see that the mean part of 
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this motion determines the position of the resonances and the oscillatory part causes 
these resonances to disappear at certain field ratios. 

Denoting these linear terms by EL{Ie,F) and on setting = e = 1, we have, to 
0{F^) 

El = -\lnIeF (^1 + ^ (23/2 + 11/^) F^ + §^(10563/6 + 772/2/2^ + 725/,t.)i^' + • • •) 

The remaining part of the Hamiltonian, FdS/dF, has a different form, eauation llBI and 
does not give rise to factors like IeF^\ this is important for the following analysis. 
Using El only in Hamilton's equations we find that <pe{t) = f^elO) — 3/„(7(t)/2 where 



g{t)= / dt 
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^ (23/2 + 11/^) F^ + §^(10563/;^ + 772/2/2, ^ 72Mt)F'> 



(24) 

The series in this integral has a finite radius of convergence, Frc{Im), so it is important 
that max(F) < Frc- We have computed this series to 0{F^'^), and used these nine 
terms to estimate Frc- By extrapolating the ratios of coefficients using Richardson's 
extrapolation we estimate that Frc ~ 0.17, 0.19 and 0.21 for Im = 0, 0.8 and 1, respec- 
tively. Using Pade approximants we obtain Frc ^ 0.18, 0.20 and 0.22, respectively. This 
provides a rough guide to the range of fields for which the following theory is valid. 

Since the field amplitude, F{t), is periodic in i, the function g{t) can be written in 
the form g{t) = gt + g{t), where g{t) is periodic with zero mean and g is the mean of El 
over a field period. With F{t) ^ Fg + F^ cos fit this becomes 

a = f, + ^(23iI + iiiI)fJf^ + If^ 



16 ^ ^^^m) ^ s 1 ^ s ' 2 

+ §^(10563/,t + 772/2/2, ^ 72blt)F, (^F^ + hF^^Fl + }^F^^]+... . (25) 

The periodic function g{t) can be expressed as the Fourier series 

p °° 

9{t) = ~^^9k sin knt 



where 



n 

fe=l 



~g, = l+^A{iul + 2ill){F^ + \Fl]+0{F' 



52 = (ll/2„ + 23/2) + 0(f4), (26) 

53 = ^i^^^(ll/^ + 23/2)+0(F4), ~gu^O{F% k>A. 

The dominant harmonic is gi. both 52 and ^3 are Q{F'^) and all higher harmonics are 
0{F'^) and may be neglected. 

Resonances in the dynamics occur when the angular frequency, 3/„g/2, resonates 
with the field frequency: the magnitude of the effect of any resonance depends upon the 
periodic component of g(t), and principally upon gi. 
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In order to see this we change to a moving reference frame, in which (j)e{t) is approx- 
imately stationary, by defining a new angle V'e = 0e + 3/„(7(i)/2 using the generating 
function F2{p,(f)e) = {<t>e + 3Ing{t)/2) p, where {ipe,p) are the new conjugate variables. 
Since g(t) is, by definition, independent of le we have p — le- In these variables the 
Hamiltonian [221 becomes 

- - EL{Ie.F{t))^ - i^/3A(J,)S(/e)sin(2V'e - 3/„5(t)) . (27) 

No further approximation has been made in deriving this Hamiltonian from defined 
in equation 1221 By definition the curly brackets contains terms quadratic and higher in 
le which are independent of tpe ; the leading term is 

E{h,F) - EUlcF) = 

and since F = F^^lsinflt we see that the terms of Km are O(F^) and 0{Fil,), and 
since fl F 0.1 (in scaled units) the mean motion of (ipeTle) is slow by comparison 
with the field oscillations: hence the relatively high frequency oscillations of £^ — El 
do not qualitatively affect the motion — a fact that has been confirmed numerically — 
and hence we may replace E — El by its mean over a field period. Retaining only the 
dominant quadratic term gives 

Km - ^ (^F! + ^F^^ Ifjl + \llF^^A{h)B{h) sin m sin (2V'e - 3/„.g(t)) . 

The second term is, for most values of 51, F^ and Eg, an oscillatory function of time 
with small mean value: in these circumstances it has little effect and may be ignored so 
Km ~ /fj giving Ie{t) ~ constant with ijjf. approximately proportional to t. However, 
for any given (il, i^^) there are particular resonant values of Eg for which the long- 
time average of the second term is proportional to sin iipe and then the nature of the 
resonant motion is qualitatively different. Near these values of Eg, Ie{t) can vary over 
a large portion of its accessible range and in some circumstances this leads to enhanced 
ionisation. 

The function g{t) is periodic and odd so we may write 

oo 

sm{2^e - 3Ing{t)) ^ ^ Jk sin {2-4;e - Vkt + kir) , Uk = 2>Ing - kVL, (28) 

k— — oo 

where the Fourier coefficients, j7/c, depend upon g^, s > 1. The coefficient J7fe is domi- 
nated by gi = 1 + 0{F^), but g2 and g^ are also 0{F^), and to this order 

oo 

Jk = Jk{zi)Jo(z2)Jo(z3) + Jo{z2)^^Js{z3) Jk-3s{zi) + {-1)" Jk+3s (zi) 

s=l 

oo 

+ MZ3)Y,MZ2) [jk-2s(zi) + {-iyJk+2s{zi) 
s=l 

CXD OO 

+ '^Js{z2)^Jr{z3) Jk~2s~3r(zi) + (- 1)"^ Jfe-2s+3r (^1 ) 
s— 1 r—1 

+ (-l)Vfe+2s-3r(zi) + (-l)'+Vfc+2s+3r(^l)l , (29) 
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where Zk = SgklnFfj^/fl. 

Using equation 1281 the mean- motion Hamiltonian becomes, 



3/'' / 1 \ 1 °° ~ 

Km = ^[F! + -F'^)l! + jlf,nF,AiQBiQ ^ J,cosi2^e - ^kt + ki:), (30) 



k— — c>o 



where the functions A{Ie) and B{Ie) are defined after equation 1221 and J} — i^Jj-i — 
Jj+i)/2, so that, to the lowest order, Jj = Jj(3F^/rio)- 

If \vj\ is small the jth term of the sum changes more slowly than all other terms, 
which may therefore be averaged over to give the resonance Hamiltonian, 

Kr = ^ [f^ + \fI^ II + \llVLF^A{h)B{h)J, cos(2^e - + jtt). 

We define the jth dynamical resonance to be at the static field, Fs^\ where Vj = 0, that 
is where coupling between Je-states is largest: the equation for Fs is 

3g{Fs, Ffj^, Im)In — or to lowest order, in scaled units, F^^^ ~ J ~ '^^ ' ' ' ■ 

We show below that near this value of Fg ionisation may be enhanced, but in section |S1 
it is shown that there is no clearly defined, precise relation between Fi"'-' and J-s'^\ the 
position the maximum in the ionisation probability seen in figures^andlSJ the two fields 
are close but the difference can be larger than the resonance width, see section in 
particuar table |21 

a) 

The position of the jth resonance, Fg ' , is, to a first approximation, independent of 
the substate quantum numbers; if one substate is ionised by this resonance others will 
be similarly affected so the effect of the resonance is not significantly changed by an 
average over substates. 

With the translation 6 = ipe — Vjt/2 + jn /2^ the resonance Hamiltonian may be 
converted to the conservative system, 



^fl - ^ [f^ + \Fi) {le - ajf + \llnF^A{h)B{h)J, cos 29 (32) 

where 



(m scaled units). 



' 3(F2 + iF2) 3(F2 + iF2) 

The Hamiltonian Kr shows that the jth resonance disappears when Jj = 0: to the 
lowest order this gives, in scaled units 

Fi^) = ijf]o and ^^0.'=) = lj-j,17o, k = l,2,--- (33) 

where j'j^. is the fcth positive zero of Jj{x). This critical value of F^^ was first derived in 
a linear quantal approximation, Galvez et al (2000), and later by Oks and Uzer (2000) 
using a Floquet approximation and by Ostrovsky and Horsdal-Pedersen (2003) using a 
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linear approximation. Recent experiments, Schultz (2003) and Shultz et al (2004), and 
comparisons with classical calculations suggest that this estimate can be inaccurate by 
up to 10%. Later we show that the present theory can be used to improve upon these 
estimates. 

The derivation of the resonance position 1311 and the resonance Hamiltonian 1321 in- 
volves a series of approximations. Before proceeding it is useful to list these in order to 
estimate the effect of ignored terms. 

1. We have used the Stark angle-action variables, defined by Hs, equation|21 For the 
action variables we use a series representation in F. For the angle variables we use 
the F = limit because these variables appear in the Hamiltonian K , eauation llSI 
only in the term which is 0{floFf^). 

2. The term dS/dF , eauation llfil has zero mean value when averaged over [ip, x) ^-^d 
because it is multiplied by the factor VIqF^j, we may approximate it by its value at 
F = 0. 

3. The mean motion Hamiltonian, Km equation 1221 is derived by averaging over 
i'n = i&i + ^2)72 which replaces sin and sin kx by Fourier series in sin 20e- The 
approximations described in points ^ and [21 truncate this series at the first term 
and approximate its coefficient to second order. 

The inclusion of higher-order terms in the mean-motion Hamiltonian 1221 introduces 
corrections 0{flFf^F) to the factor nFf^A(Ie)B{Ie) and adds further terms corresponding 
to the harmonics shi2p(j)e, p = 2, 3, • • •. Crucially this means that the estimate, Vj = 0, 
of the j resonance position is not affected by the approximations made: that is, the 
position of the dynamical resonance is determined solely by the parts of the Stark 
Hamiltonian, E{T), equation 1101 which are linear in /g. 

A better estimate of Fs''^ is therefore obtained using the series [23 which has been 
evaluated, using computer assisted algebra, to 0(F^^). In section we use this to 
obtain better estimates of Fs and the values of F^ at which the resonances disappear. 

4 Resonant lonisation 
4.1 Qualitative discussion 

Here we show how the dynamical resonance described in the previous section can en- 
hance ionisation. The connection is qualitative, but explains many featues of the loni- 
sation probablity. 

The Hamiltonian Kn, equation 1321 is similar to that of a vertical pendulum, but 
there are two significant differences. First, le is confined to the region |/e| < /„ — /,„ 
with natural boundaries at 1^ — ±(/„ — Im), where B{Ie) = 0, see equationl22l Second, 
the coefficient of cos 26* depends upon /g. The fixed points of Kji are at the roots 
of dKfi/dle = dKu/dO = 0, and for this analysis it is convenient to replace J'j by 
|j7j |; when J} < this represents a physically unimportant translation in 9. The 
equation dKn/dO = gives, for < 6* < tt, 6* = and 7r/2. At = 7r/2 the equation 
dKn/dle = has a single root near 1^ ~ ctj and this fixed point is a centre. At = 
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there are generally three roots: a saddle near /g = Uj, but there are two others near 
le = ±(/n — /rn), becausc of the square-root singularity in A{Ie)- If the phase 
plane is projected onto a sphere with latitude -0, so le — {In — Im) cosV', it is seen that 
there are phase curves with centres close to, but not at, the poles and which enclose 
the poles: in the Cartesian coordinate system [9, le) this produces the two extra fixed 
points. The physically significant fixed points are near /e = aj and these exist only if 
\aj\ < In — Im, approximately. 

By plotting the contours of Kb. for fixed and and with Fs increasing so aj 
increases from below — (J„ — Im) to above /„ — Im we see how the resonance develops 
and why, in certain circumstances, ionisation is enhanced if Fs ~ Fs ■ 

The following five figures show the contours of Kr, near the j = I resonance for 
Im = 0.2, SIq = 0.0528 and — 0.13, corresponding to figure |21 For these graphs 
we use Kfi with gi and g given by equations 1261 and 1251 fto 0{F^)) respectively; these 
give vi ^ when i^i^'' = 0.0168. In each figure Fs and F)j are fixed so, according 
to the adiabatic ionisation criterion, there is a critical value of le, above which orbits 
ionise, given by the solution of Fcrit(/e,^m) — Fs + F^: for the parameters used here 
the critical value of changes from le = 0.50 {Fs = 0.0152) to le = 0.42 (F, = 0.0183). 
The maximum of Fcrit{Ie,Im) is at /e = 1 - so if F^ -|-F^ < Fcrit(l -\Im\,Im) there 
is no ionisation, even at a resonance. The upper solid horizontal line in each figure is at 
the value of le at which Faitile, 0.2) = Fs + Ff^, so orbits straying above this line will 
ionise. The lower solid line is le = —0.4 is taken, for illustrative purposes, to be the 
initial state. Note that in this case F^ — Fs < 0.13, so the other adiabatic boundary 
defined in equation 12 II does not lead to ionisation. 




0.5 e/7t 1 0.5 e/jt 1 0.5 e/n i 

Figure 6 Figure 7 Figure 8 

ai = -0.75, Fs = 0.0152 «i = -0.25, Fs = 0.0163 ai =0, Fs = 0.0168 




0.5 e/7I 1 0.5 e/TT 1 

Figure 9 Figure 10 

ai = 0.25, Fs = 0.0173 ai = 0.75, Fs = 0.0183 



As Fg increases from 0.0152 to 0.0183, through fJ^'' = 0.0168, the centre of the resonance 
island moves upwards. The physical effect of this is understood by considering a field 
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suddenly switched on, with the initial value of <j)e uniform in (0, 27r) and the initial value 
of /e, to be -0.04. 



• FigureEl Fs — 0.0152 {ai = —0.75). The adiabatic condition shows that orbits 
for which /e(i) > 0-50 ionise; at this field there is no resonance island and no 
ionisation from the initial state. 

• Figure Ul Fg — 0.0163 (ai = —0.25). The resonance island exits; it intersects the 
initial state, but does not overlap the ionising region, so there is no ionisation. In 
practice the demarcation between ionising and non-ionising regions is less sharp 
because the averaging approximations used to derive this simple picture replaces 
unstable manifolds by separatrixes. 

• Figure IHI Fs — Fi^'' — 0.0168 (ai = 0). The resonance island is in the centre of 
the phase space. The ionisation criterion is practically the same as in figure El so 
orbits with leit) > 0.50 ionise. Now, however, the resonance island can transport 
initial states to the ionising region, /g > 0.46. Note that not all orbits trapped in 
the resonance island will ionise, but only those near the separatrix. We shall see 
in section |S1 how this affects the ionisation times. 

• Figure Fs = 0.0173 (ai = 0.25). The centre of the island is now at 4 = 0.2 
and its separatrix just dips below the initial state, so few orbits ionise. In these 
circumstances it is shown in section|Slthat the ionisation probability can be affected 
significantly by the way the field is switched on. 

• Figure mi Fs — 0.0183 (ai = 0.75). As for figure El the island no longer exists 
and there is no ionisation from the initial state. 

This qualitative description of the ionisation process suggests that for a microcanonical 
distribution of substates the background ionisation increases as Fs increases across a 
resonance, because decreases, as shown in figure ITU 

The centre of the jth resonance island is at approximately 



In this field range a high proportion of initial values of /e(0) may lead to ionisation: 
outside this interval the resonance does not exist. 

This qualitative description of the ionisation mechanism shows that for a system 
initially in a given /m-substate there are several conditions necessary for a resonance to 
enhance the ionisation probability. 

1) The field amplitudes must be sufhciently large that there is ionisation for some 
value of /e, for the given I^- 




(in scaled units). 



see equation 1321 so it exists only for Fg in the interval 




■f^l + ^Fl 



(34) 
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2) The field amplitudes must not be so large tliat Pi = 1. 

3) If Fci.it(/e , Im) = Fs + Ffi, then the island width must exceed — /e(0), otherwise 
the initial state cannot be transported to an ionising state. There is, of course, 
a similar relation for the boundary defined by Fcrit(— ~Im) — Ffi — Fg when 
i^M > Fs. 

The first two of these conditions define a region in the {Fs , F^)-plaTie in which the 
resonance may enhance ionisation. The boundaries of this region depend upon and 
are the complement of the region defined by the two conditions Pi = and Pi = 1 for 
all le- Using the adiabatic assumption, Pi — 1 if Fg + F^^ > max(Fcrit) and Pi = if 
Fs + Ffj^ < min(Pcrit) and F^ — Fg < min(Fcrit), if Fg < F^, see equation!^ In the case 
Im — 0.1 this region is shown by the shaded area in the following figure: outside this 
region the resonance can have no effect. If > 0.1 the equivalent region lies inside 
the area shown. 



0.25 




0.05 0.1 0.15 0.2 0.25 



Figure 11 Diagram showing the regions where Pi = and Pi = I 
for Im = 0.1. In the shaded area < Pi < 1: only in this region are 
the resonances potentially visible. 



Inside the shaded region a resonance affects the ionisation probability only if condi- 
tion 3 above is satisfied. Below, but near, the upper boundary P,; ~ 1, so resonance peaks 
are barely noticeable, see for instance the j = 5 resonance in figureHat Fg ~ 0.148. The 
break down of adiabatic invariance broadens these boundaries but in a manner difficult 
to estimate, although the effect increases as fio increases; tunnelling also affects these 
boundaries. 

There are three parameters of the resonant island that affect Pi{Fg). These are most 
easily estimated by setting Fg = Fs-'\ so aj = 0, and /,„ = as well as using the 
lowest-order estimates of all variables: these values are used in the remainder of this 
section. 

The first parameter is the island area, Aj, which determines how the classical reso- 
nance affects the quantum dynamics. An estimate for this is, 



A- 2/„ . ^1 snoF^lJ', 

— ■ sm ' 



1 z?2^ 1 ACi IT I PI' "^1 ~ "^jU^f^/^s), Fg — F^^h 



2TT TT y 3(F| + iF2)+4f^^J^^|J/| 

the form of the ] 

is proportional to y/lJ'jl- For the parameters of figure^ ^Iq — 0.098, = 0.1 this 



For other values of and aj the form of the resonance Hamiltonian shows that Aj 
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gives Ajl2T: = (0.80, 0.05, 0.35, 0.25)/„ for j = 1, 2, 3 and 4. In this case no = 43 so 
the approximate number of states associated with these islands are 34, 2, 15 and 11, 
respectively. Thus if no is decreased, all other scaled variables remaining the same, we 
should expect, in quantum dynamics, the j — 2 resonance to become less prominent 
than the other resonances. 

The second parameter is the island width, A/e, that is the maximum distance be- 
tween the two branches of the separatrix. It is difficult to derive a simple estimate of 
A/g; here we simply note that it is proprtional to | Jj|, defined above. A necessary con- 
dition for enhanced ionisation is that A/g is larger that — /e(0); otherwise transport 
to ionising regions docs not occur. Notice that this condition is independent of the 
principal quantum number, n, unlike that discussed above. 

The third important classical parameter is the period of the mean motion inside the 
island; as we shall see, this determines how rapidly a resonance develops, section|5l and 
how it is affected by the field envelope, section The frequency, tUj, of the motion 
inside the resonance island is approximated by expanding the resonance Hamiltonian, 
equation 1321 about its centre and near the island centre we obtain, in scaled units 

Fs = ^jno. (35) 

This estimate gives the largest value of the frequency in the island; for motion nearer 
the separatrix cuj is smaller. 

In this section we have shown how a dynamical resonance can enhance the ionisation 
probability and have derived some approximate necessary conditions. The analysis uses 
the resonance Hamiltonian, Kr, equation 1321 derived using two stages of averaging. 
Moreover, an important part of this Hamiltonian is the factor J7j, which, for small Qq, 
oscillates between its maximum and minimum values for relatively small changes in . 
Hence, whilst Kfj provides a good qualitative description, for any fixed {^OtF^) the 
details may be wrong; for instance the field at which a resonance disappears is given 
inaccurately by this approximation if f^o is small. 

4.2 Resonance positions 

In this section we examine ionisation from a particular substate and compare theo- 
retical predictions with exact numerical calculations. We choose the low frequency 
ilo — 0.011414, (no = 21) (to minimise non-adiabatic effects), fix /,„ = 0.2, use the 
initial condition Jg = —0.4 (so there is no average over substates) and put = 0.13. 

Since le is an approximate constant of the motion and Fcrit(— 0.4, 0.2) = 0.1984, 
if the dynamics were adiabatic we should expect complete ionisation when exceeds 
Fcrit — Ff^ = 0.0684 and no ionisation for smaller static fields. At a resonance Ie{t) varies 
over part of its accessible range and since min(Fcrit) = 0.1357 we might see the effect 
of the jth resonance if Fs > 0.0057, provided the size of the the resonance island is 
sufficiently large. This simple analysis suggests that the resonances 2 < j < 18 could 
be seen via ionisation: in practice, for reasons to be discussed later the j = 2, • • • , 6, 11 
and 14 resonances are not observed. 

In figure El we show the classical ionisation probabilities for the envelope 16-50-16 
in which the j =7-10, 12, 13 and 15-19 resonances are clearly visible, but the j = 6, 



Ft 
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11 and 14, marked by the arrows, are missing: other calculations show that the j = 5 
resonance is also missing and theory suggests that the j = 2-4 resonance islands are too 
narrow to affect the ionisation probability, that is Ale < le ~ -^e(O), as discussed in the 
previous section. 




0.02 0.03 0.04 0.05 0.06 0.07 0.08 

Figure 12 Ionisation curve for Qq = 0.011414, (no = 21), = 
0.13, with initial conditions I2 = 0.2, Im = 0.2. 



In table 12 are listed some parameters associated with the j — 7-15 resonances. Here 
the resonance width, AFs''\ is defined to be the difference F+ — F~ where are 
respectively the smallest and largest values of Fs on either side of Ts at which i-^ = 0. 
Notice that this width is generally less than the difference \Fs''^ — J-s''^\. In these calcu- 
lations 500 orbits were used. The values of J-'i''\ the static field at which Pi is largest, 
are computed using a grid 10~ in Fs, and J-s is taken to be at the maximum value 
of Pi on this grid. The value of Fs''^ is the root ofg{Fs, Ff^) = jfto/S: below the double 

lines Fs + Ffj_ > Frc, the approximate radius of convergence of the series for see 
equation!^ 

The function g(Fs,F^) is known only via its series expansion, equation 1251 which has 
been computed to 0{F^'^), that is the first nine terms. From the discussion after equa- 
tion |^ since max(i^) = Fs + F^, we expect any theory based on the series represen- 
tation of g to be valid for those resonances satisfying F^^^ < 0.18 (/,„ = 0) and 
0.22 (/,„ = 1). 

It transpires that if F is near the upper boundary the values of Fi"'-* are sensitive 
to the number of terms in the series for g and extrapolation is necessary to estimate 
converged values. Here we consider two methods of extrapolating and give reasons 
which suggest that the Pade approximant is more reliable. All the following results 
are obtained by substituting F^ = 0.13 and Im = 0.2 into the series for g and then 
manipulating the resultant power series in F^: for completeness we give this series: 

= 1.057 + 0.08759a: -f 0.1079x2 -I- 0.1073x3 -f 0.05903a;'^-f 0.01588a;'^ 

+0.001861a;^ + 7.700 x lO^^x^ + 6.509 x 10" V, x = (lOFs)^. (36) 

(7) 

For the j = 7 resonance the lowest order approximation gives Fg ~ 7V,o/3 — 0.02663, 
which is about 7% too large. Eight other estimates can be obtained by truncating 
the series for g{Fs)/Fs at F"^^ , k = 1, 2,- ••,8: these are 0.02566, 0.02536, 0.02522, 
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Table 2: Table showing the positions, heights and widths of the resonances seen in figures [T^ 
For completeness, the values of F^^^ are 0.00360, 0.00720, 0.01079, 0.01438, 0.01795 and 0.2151, 
for j — 1-6, respectively. 
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0.02841 
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0.03544 
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0.04235 


0.04240 


0.73 


32 X 10 
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13 


0.04581 


0.04576 


0.73 
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-5 


14 




0.04906 





< 10-^ 




15 


0.05269 


0.05229 


0.94 


77 X 10 


— 5 



0.02515, 0.02511, 0.02509, 0.02507 and 0.2506. This sequence appears to be converging, 
but has not reached its limit (to the four significant figures quoted). Suppose we have 

M estimates Tp for Fs''\ using p = 1, 2, • • • , M terms of the series OSl then using the 
method of Richardson we assume that Tp = Fi^'^ + E^li^ Arp-''. These M equations 
may be solved for the unknown Fi"'"' to give the estimate fP ~ 0.02505, which differs 
from Ts''\~ 0.02495) by 10"''. Despite this relatively small inaccuracy we note that a 
Monte-Carlo calculation (with 500 orbits) gives Pi(0. 02505) = 0, with P,{Fs) ^ in the 
interval 0.02495 ± 5 x 10^^, which does not overlap with our estimate of Fs'^\ 

Another approach is to form a Fade approximant oig{Fs)/Fs using the expression l36l 
treated as an eighth degree polynomial in F^. The coefficients, and hence the positions 
of the poles, of these approximants depend upon and F^: we find that the position 
of the pole nearest the origin is relatively insensitive to /„, but changes significantly 
with Ff^. Therefore for < 0.09 we use a [2/2] approximant (in F^) and for > 0.09 
we use a [3/3] approximant. 

For the case considered here, F^ — 0.13, Im — 0.2, the relevant Fade is 

g 1.057 - 9.542a; + 64.89x2-156.4x3 .„^^2 

— = K 5—, X = (lOF,)^. (37) 

F, 1 - 9.858a; + 59.36a;2 - 198.0x3 ' ^ ' ^ ' 

With this approximation for g(Fs) the equation Sg — jQo gives the values of Fi"''' quoted 
in table El 

(15) 

For the j = 15 resonance the zero-order approximation to Fs is 0.0571, which is 
about 8% too large. The eight other estimates, obtained using the truncated series for 
giFs), equation EHl are 0.05479, 0.05395, 0.05347, 0.05313, 0.05287, 0.05265, 0.05246 
and 0.05229. Richardson's extrapolation gives Fi^^^ = 0.04858, which is 8% smaller 
than J^s''\= 0.05265), and no improvement on the zero-order approximation. The [3/3] 
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Pade approximant, cauation l37l gives Fs — 0.05229, approximately 0.7% smaller than 
In this case the Pade approximant seems to provide a more reliable method of 
extrapolating the truncated series for 'g{Fs). We note that this resonance field is on the 

(15) 

edge of the validity of the series expansion, so any estimate of Fs based on the series 
may not be accurate; in these circumstances, however, the Pade approximant is more 
likely to provide an accurate estimate of the exact function. In all future estimates of 
we therefore use the Pade approximant. 
We now turn our attention to the resonances missing from figure [T^ For this analysis 
we use the resonance Hamiltonian used to plot the contours in figures HillOl which, for 
the reasons discussed at the end of the previous section, provides only a qualitative 
description. 

Using the simple approximations for g and gi, eauations l25l and l26l we find than the 
width of the resonance island is proportional to 

A{3)^j;(^^~gi{F^'\F,)^ where g{F^^\ F,) = ^jU^. (38) 

An overview of the widths of the j — 2-15 resonances is given by the graph of this 
function with j taking all real values. This graph is shown next, with integer values of 
j being marked by the circles. 
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Figure 13 Graph of the function A{j) for = 0.13, Vlo = 0.011414 and 
/m = 0.2 



We relate this graph to the ionisation curve in figure [T^ bv recalling that a dynamical 
resonance affects the ionisation probability only if it can transport an orbit to a region 
le > lei where is defined by Fertile, Im) — F^^ + Fs''\ so the resonance island width 
must exceed the difference — le (0) . 

The values of ^i^i"''^ decrease with increasing j; for j = 2, — 0.75 and for 

j ^ 5 — 0.43, but for these cases the maximum possible size of the resonance island is 
smaller than — /e(0). For j = 6, 11 and 14 the resonance island is seen from figure [T^ 
to be very small so these resonances do not affect the ionisation probability. 

For j = 7-10 the simple Hamiltonian with contours shown in figures IfillOl suggests 
that the resonance island is slightly too small for enhanced ionisation. But non-adiabatic 
effects, the approximations used — recall that two averaging approximations have been 
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used to derive this Hamiltonian — and the field envelope will broaden these boundaries. 
The J = 11 and 14 resonance islands are predicted to be too small to promote ionisation, 
whereas this simple approximation predicts enhanced ionisation for all other j. 

4.3 Resonance disappearance 

The jth resonance has no effect on the dynamics at those values of where J'j — 0, see 
equation 1321 With g(Fs,_F'^) approximated by a Pade approximant, as in equation 1361 
the equation g = jTJo/3 provides an expression for Fs''\Ff^) and then equation EHl can 
be used to obtain the numerical value of , associated with jj^,. In the following 

table we give values of Fs and F//' , nearest 0.13, at the j — 5-15 resonances for 
riQ = 0.011414. In this example the difference between and f//'''^ is about 0.01. 

Table 3: Table showing the value of _F'i"'' and F^f'^^ nearest 0.13, at which each resonance 
shown in figure WI\ where -F^ = 0.13, disappears. 
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10 


pi,) 


0.0215 
0.130 


0.0249 
0.135 


0.0287 
0.128 


0.0319 
0.133 


0.0358 
0.126 


j 


11 


12 


13 


14 


15 


pU) 


0.0390 
0.130 


0.0420 
0.135 


0.0460 
0.127 


0.0489 
0.131 


0.0515 
0.136 



Observe that for j = 6, 11 and 14, F^f' ' ~ 0.13 and that for j = 8 it is close to 0.13 
and at this resonance max(Pi) is relatively small. 

In the experimental results reported by Galvez et al (2000) it was shown that res- 
onances disappear at certain field values {J-i''\ F^'''''^), given approximately by equa- 
tion 1221 Since then two sets of more accurate measurements have been made. First, 
disappearances in the 8.105 Ghz cavity, with scaled frequencies mostly in the range 
0.0731 < r^o < 0.136, are reported by Galvez et al (2004). Second, Schlultz (2003) has 
reported results for a cavity with frequency 3.5539 GHz, and scaled frequencies in the 
range 0.035 < ^Iq < 0.16, see also Schultz et al (2004). In those papers experimen- 
tal values of {!Fs^\ F^i''^'^) are compared with classical Monte Carlo estimates: here we 
compare some of these with theoretical values. 

With the Monte-Carlo method it is not feasible to compute the exact values of 
{Fy \ F^"'"'"''), so we determine an interval for each of Fg and F^ in which the resonance 
cannot be distinguished from the statistical fluctuations. The mid point of this rectangle 
is taken to be the point of disappearance. 

For the 29 cases considered by Schultz (2003) the relative difference between the 
experimental and computed values of Ts''^ is less than 1% for 22 cases and between 1% 
and 2% in 5 cases: for Fli'''^^ the corresponding percentages are 11 and 16 respectively. 
For a comparison with theory we compute {Fs^\ F^^'^'') for between and 1 and 
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compare the average with the experimental values of {Ts , F^''' ). Of the 32 cases 

considered the relative difference between Fs''^ and J^s''^ is less than 1% for 18 cases and 
between 1% and 2% in 10 cases: for Fjf'''^ the corresponding percentages are 24 and 6 
respectively. 

Some typical comparisons between the Monte Carlo calculations and the theoretical 
values, for various scaled frequencies, are shown in table 01 



Table 4: Comaprison of the theoretical values of {Fp , Pji''''^) and the Monte Carlo estimates 
of (Ti'\F^'-''-^) for various scaled frequencies and values of j. In the extreme right column is 
J'^P' +F^'''^^ because the discussion after eauation l24l suggests that this theory will be unreliable 
if this value exceeds 0.19. 
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4.4 Variation of resonance position with 

Now consider the variation of Fs''^ with /,„. From the series for equation we 

see that Fs^^ decreases as /„i increases, but that the difference between the largest and 
smallest values is approximately lljTio(9-F'^ + (jf2o)^)/432: for the data in table|21this 
gives 4.4j X 10"^ approximately. The results obtained using the Fade approximant for 

g are given in table [S| Columns 2 to 4 show the mean (^Fs''^'^, averaged over /,„, and 

the minimum and maximum values of Fs , for Ff^ — 0.13 and f2o — 0.0114. 
This data shows that, for this low frequency, the above simple estimate of the spread is 
reasonable for j < 10, and that it is comparable with the difference between the Ts 
and F^^\ given in columns 2 and 3 of table |21 and to the resonance width, column 5. 
Note also that as j increases, and Fj-'-* + F^ tends towards the radius of convergence of 
the series for the difference decreases: it is not known if this effect is real or due to 
the approximations used. 

Although the position of the dynamical resonance is very weakly dependent upon 
Im its effect on Pi{Fs) can depend strongly upon /„, because — that is the solution 
of Fcntm,Im) — Ff^ + Fs''^ — depends upon /,„: for instance with j = 10, varies 
from 0.04 (/„ = 0) to 0.3 (/„, = 1). 
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Table 5: Table of the theoretical resonance positions for the parameters defined in table |51 
Ff_i — 0.13 and f^o = 0.0114. In column 5 is the difference between the largest and smallest 
values of F]^'. 
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4.5 Resonance widths 

The shapes of the resonances seen in the experimental data and the classical simulations 
are complicated. In particular the resonances are not normally symmetrical about J-g , 
with details depending upon Ff^ and j. Here we show show the classical widths can be 
defined and computed and discuss some of the reason for the shapes observed. A detailed 
comparison between the classical simulations and the experimental data is provided in 
the companion paper, Galvez et al (2004), and there it is shown that the experimental 
results display the same complexities. 

In order to estimate the widths it is necessary to isolate the resonances from the 
background. This is achieved by noting that on either side of the resonance Pi{Fg) 
increases approximately linearly, figure rPTl below. The background may be eliminated 
by subtracting these straight line segments from Pi(Fs) to give an adjusted probability 
that is approximately zero on both sides of the resonance, as shown on the right of 
figure [CT The only complication with this procedure is that the straight line segments 
have different gradients, so we form a new fit to the background, = m{Fs) F^ + c{Fs) , 
where the gradient, m{Fs), and constant, c{Fs), change smoothly between the values 
either side of the resonance, (rni,ci) and (^2,02) respectively. If the straight line 
segments are on the intervals Fi < Fg < F2 and F3 < Fg < F4, with F2 < F3, (chosen 
by eye) we set 

m{F) = — - + —-^ — ^tanh(aF + /3), a= , (3- aj 



F3 — F2 F3 ~ F2 

with a similar fit for c{F). 

In the following two figures we show how this process works for the j = 2 resonance 
with no = 47 (Oq = 0.1278) and F^ — 0.1, which is a fairly typical example: in these 
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calculations a microcanonical distribution of initial states with 1296 orbits, for each 
value of Fs, is used and the envelope is 16-113-16. 



J I , , , , , , , , , , , , , , , 1 Q2 




U.3 i ' ' ' 1 ' ' ' 1 ' ' ' 1 ' ' ' h I 1 ' 1 • 1 ' 1 ' 1 1 

0.06 0.07 0.08 0.09 0.1 0.076 0.078 0.08 0.082 0.084 

Figure 14 A graph of the j = 2 resonance for Qq = 0.1278 (no = 47) with = 0.1. On the 
left is the ionisation curve with the straight hne fits, as described in the text. On the right is the 
dilTerence, having subtracted the background. 



The left hand panel of figure H"!] shows that either side of a resonance Pi{Fs) is ap- 
proximately linear, but with different gradients. The difference between these gradients 
changes with F^^ and, of course, is zero when the resonance disappears. The right hand 
panel shows the graph of Pi{Fs) — P^^{Fs) which highlights the resonance shape; this 
is clearly asymmetrical about the maximum. The graph shown is typical though the 
degree of asymmetry changes with j and F^. The position of the resonance, the two gra- 
dients mi and TO2 and the width of the adjusted ionisation probability all provide tests 
for any theory. Comparisons of these parameters obtained from classical calculations 
and experiment are given in Galvez et al (2004). 

(2) 

In this particular example JFs = 0.0802: the width at half-height is about AF, = 
0.0018, with the base being nearly 4 times wider, 0.007. The calculated position of this 
dynamical resonance varies from Fs — 0.0797 (/„ — 1) to 0.0800 (/,„ = 0). The spread 
in Fi^-* due to the variations in is therefore about 0.0003, about ^ of the ^-height 
width seen in figure FPTI this difference is fairly typical. 

We now consider some of the factors determining the resonance shapes, and show how 
these are partly determined by a combination of substate averaging and non-adiabatic 
effects. Consider ionisation from a given substate: for illustrative purposes choose 
no = 47 {flo — 0.1278), /,„ = 0, = 0.1 and use a 4-50-4 envelope with initial 
conditions /e(0) = —0.9, —0.8, —0.6 and —0.4. The values of Fdit — F^ are depicted by 
the arrows in each of the four graphs and adiabatic invariance suggests that = 1 to 
the right of these arrows. 
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0.05 0.1 0.15 0.2 0.05 0.1 0.15 0.2 

Figure 15 lonisation probabilities for Qq = 0.1278, Ffj, = 0.1, Im = and /e(0) = —0.9 
and -0.8. 



When /e(0) = —0.9 the j — 1-5 resonances are clearly present; the positions and full- 
widths, as defined in section are given in tableEl For 1^ ~ —0.8 only the first three 
resonances are visible because Fcrit = 0.256. In the next two figures we see the j = 2 
resonance disappearing as /e(0) increases. 



I_^(0)=-0.6 I^(0)=-0.4 




0.02 0.04 0.06 0.08 0.1 0.12 0.02 0.04 0.06 0.08 0.1 0.12 

Figure 16 lonisation probabilities for Qq = 0.1278, = 0.1, Im = and /e(0) = —0.6 
and -0.4. 



Table 6: Values of F^', T^' and the full resonance width, for the resonances shown in 
figures \TE\ and \Wi the latter two items being computed as in table |5| with the fs-grid being 
10~^. The radius of convergence of the series 1^ for this problem, is about 0.18. 



j 


12 3 4 


1 2 3 




le = -0.9 le = -0.8 


pin 

full- width 


0.0413 0.0801 0.109 0.122 
0.0407 0.0804 0.119 0.154 
0.0014 0.0018 0.0028 


0.0413 0.0801 0.109 
0.0407 0.0805 0.1185 
0.0017 0.0027 0.009 







le = -0.6 


le = -0.4 






0.0408 


0.0804 


0.0408 0.0788 




full-width 


0.0026 




0.0035 
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In this example adiabatic invariance predicts that the j = 4, 3, 2 and 1 resonances 
disappear when 1^ — —0.80, —0.60, —0.20 and 0.63, respectively: figures [TslandlTBlshow 
these predictions to be approximately true. This data also shows that as le increases 
and the adiabatic boundary encroaches upon each resonance it broadens, acquires an 
asymmetry and eventually disappears. In other cases, when F/j, > Fg, an ionisation 
boundary can also encroach from the left, figure[51 and this will also change the resonance 
shape. 

This example, which is typical, suggests that the width and shape of the micro- 
canonical averaged resonances is caused mainly by the effect of the separatrix between 
bound and free motion, which distorts nearby resonances, rather than the variation in 
the resonance position with /,„ . It is therefore difficult to provide theoretical estimates 
of the resonance width and shape. 

5 Resonance time-scales 

The classical adiabatic ionisation mechanism, described in section ISTTl suggests that, in 
the absence of a resonance, ionisation occurs when F ~ X{t){Fs + Ff^) reaches a critical 
value defined by the condition F — Fertile , Im)'- at this time ionisation from a particular 
orbit occurs within a Kepler period. This behaviour has been checked numerically when 
i^s = (Rath, 1990) and the results of the present calculations, where F^ > 0, show the 
same behaviour. With increasing scaled frequency the dynamics becomes less adiabatic 
although this behaviour persists, albeit with the boundaries becoming blurred; this 
behaviour is seen clearly in figure 2 of Richards (1996a). 

The classical ionisation mechanism at a resonance is different because here for an 
orbit to be ionised it must first be transported into a region of large le , that is smaller 
-fcrit , by motion round the resonance island. Hence the rate of ionisation will depend 
upon the period of this motion; at the island centre this is given approximately by 
equation which shows the period to be 0{F~^/'^^q^^'^). Thus we should expect the 
time dependence of the ionisation probability on and off resonance to be quite different. 

These predictions can be checked by computing the time at which Pi {Fg , t) reaches a 
given proportion of its final value. In the following figure we show two graphs which allow 
comparison of this ionisation time with the ionisation probability. The upper graph is 
the ionisation probability, Pi{Fs), for fig = 0.0528 (no = 35), = 0.13 starting in the 
initial state (Jg,/™) = (—0.4,0.2), using a 16-50-16 envelope and 1600 orbits, which is 
the same as in figure |21 the lower graph is the time, T/j, at which Pi{Fs,t) reaches half 
its final value: with Th is measured in units of the field period: the horizontal line is at 
T = 16Tf, the time when the field amplitude reaches its maximum. 
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Figure 17 lonisation probabilities, upper graph, and ionisation times, lower 
graph, for the parameters defined in the text. Here Qq = 0.0528 and = 0.13. 



There are several features of this comparison worthy of note. 

(i) For Fs > 0.06 ionisation occurs close to end of the switch-on time, T = 16Tf. 
Since Fcrit{— 0.4, 0.2) — 0.198 the adiabatic condition suggests, that away from 
resonances and these initial conditions, bound states exist only for < 0.07. 

(ii) At the j — 1, 2 and 3 resonances, ionisation occurs some time after the field has 
reached its maximum amplitude, with the longest delay occurring at the edges of 
the resonance and the shortest near the maximum in Pj. This is consistent with 
the description given in section^ where it is shown that close to the resonance 
edge transport is near the separatrix where motion is slowest. The formulaic for 
ujj, gives, for these parameters, the 1/2-period near the resonance island centre of 
about 14T/, which is consistent with the lower graph of figure ITTI 

(iii) The local maximum in Tfi at Fg ^ J^s cannot, at present, be explained. 

(iv) The ionisation time near the local maximum in Pi at F^ = 0.0428 has the same 
shape as those near the j — 1 and 2 resonances, but ionisation clearly takes 
longer, suggesting that this structure is due to a higher-order resonance. Linear 
interpolation between the j = 1-3 resonances suggests this could be the j = 2| 
resonance. A similar calculation suggests that local maximum in at Fg ~ 0.0537 
could be the j = 5-^ resonance. 

It is not easy to see what produces these non-integer resonances. Second-order 
perturbation theory applied to the mean motion Hamiltonian 1221 does not appear 
to give 1/3 resonances; this suggests that higher harmonics of (pe are required and 
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these occur, at this level of averaging, only if higher order terms of diiipjX) s-^d 
6*2(1/', x) I equations [T^ and 1^ are included. 

(v) For Fg > 0.05 the boundary at i^crit— — 0.068 is begining to affect the dynamics 
and seems to be interfering with the j = 3 resonance. 

6 Envelope effects 

In the previous section it was shown how resonance islands affect ionisation times. Here 
we examine the effect of the envelope switch-on time, — on a particular 

resonance. At this point it is useful to recall that the dynamical resonances discussed 
here are unusual because each exist only for a narrow range of Fg and within this 
interval the resonance island moves from the lower to the upper edge of phase space, see 
figures I6I10I Moreover, the motion inside a resonance island is very slow, equation 1351 
If Fg ~ i^i"'' , for some j, then for most of the switch-on period the resonance island does 
not exist. But for some time close to the island develops at the bottom edge of phase 
space and as f — s- it moves up through phase space and through the initial phase line. 
As this happens the line is distorted, with the amount of distortion depending upon the 
relative values of dX/dt and the frequency of the motion in the island. For short switch 
times the initial phase line /e(0) evolves into a nearby line at t = Ta, as shown in the 
left panel of figure 1201 below. For relatively long switch times there may be enough time 
for the initial phase line to develop an incipient homoclinic tangle and become quite 
complicated, as seen in the right panel of figure 1201 The examples considered next show 
how changes in the switch-time, Ta, can dramatically affect the ionisation probability. 

The demonstration of this effect is in two parts. First we show some exact numerical 
results illustrating how the ionisation probability changes with Ta- In this example we 
chose the low frequency, Hq — 0.011414, (ng = 21), = 0.2, F^ = 0.13 and examine 
Pi{Fs) in the vicinity of the j = 7 resonance for the envelope Na-50-Na, for Na = 1-40. 
From table[21we see that when Na = 16, Pi{Fs) has its maximum at = 0.024950 and 
that the resonance width is AFs ~ 8 x 10~^. 

In the following two figures we show how Pi {Fg ) changes with Na : for these calcula- 
tions we used 900 orbits and a grid SFs — 2 x 10^®, so there are 50 data points for each 
unit of / = {Fs - 0.024925)101 

The first figure shows ionisation probabilities for Na — 1-17. The variation of Pi 
with Na shows a surprising amount of variation; in particular we note that for Na — 15 
the ionisation probability is zero across the resonance. 

Also observe that JFp-' changes by AJ^ ^ 5.2 x 10^^ for Na in this range, and that 
this is comparable to the resonance width. This explains why an unambiguous relation 
between Fs and Ts does not exist. 
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Figure 18 Some graphs of the ionisation probability in terms of the scaled field 
f = (Fs - 0.024925)10* across the j = 7 resonance for = 0.011414 and switch times 
for 1 < Afa < 17. 



In the next two figures are shown the ionisation probabilities for 23 < Na < 40. With 
these longer switch times more structure is seen. For instance, with Na — 31 and 36, 
Pi{Fs) has two local maxima and for both Na = 31 and 37 the probability has a long, 
low plateau after the maximum. A qualitative explanation of these feature is given next. 



0.4 



0.2 



P.(F) 


/-M N=24 




/ ^27 \ 




J-< ■ \ 




















if 11 n;. V 



-0.5 



1.5 



0.6 


P.(F) 




\n=40 


0.5 








0.4 




31,, /3I/ - 


--,36 \ 


0.3 














.,.32 \ \ 


0.2 


36 






0.1 




X^v'' 37\ \ 













-0.5 



0.5 



/ 



1.5 



Figure 19 Some graphs of the ionisation probability in terms of the scaled field 
f = (Fs - 0.024925)10* across the j = 7 resonance for Qo = 0.011414 and switch times 
for 23 < A^a < 40. 



The behaviour depicted in figures ^1 and ^| can be understood qualitatively using 
a combination of the mean- motion Hamiltonian, equation 1391 below, and the resonance 
Hamiltonian. We assume that initially the system is in a given /g-state with its conjugate 
variable uniformly distributed in (0,7r). This initial phase line, Co, evolves during the 
switch-on period in the mean- motion Hamiltonian: in scaled units, with = 1, this is 

^ -3AW^/ 1 \ 1 
16 I ^ 2 ^/ 4 



^™ = + T^Fu I; + -\{tW^A{I,)B{Q sin (2^^, - 3g(i)) sinf^o*, (39) 



where g(t) = J^dtX{t) {Fg + Ff^ cos ilot) . For t > Ta (and before the switch-off time), 
A = 1 and this Hamiltonian simplifies to 

3 / 1 \ 1 / 3F \ 

K,n = (f^ + j ll + -^^oF^A{QB{h) sin f 2V', - 35, - iF,t + sm^ot\ sinQnt 
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where ga — g{Ta) —TaFg ^ sin f^oTa- Near the jth resonance this can be approximated 
by the resonance Hamiltonian, 



1 /3F 

-r!oF^A(/e)B(/e) Jj i-^ ] cos{2e, + jtt), (40) 



where Uj is defined after equation 1^ and 29e — 2ipe — {3Fs — j^o)t — Sga, {t > Ta). 

During the period < t < Ta the initial phase curve evolves according to the 
Hamiltonian KWI into the line Ca- Hence by plotting the line Ca and the contours of Kj 
(the dashed lines), we obtain a qualitative picture showing how Pi can be affected by 
the field switch. 

In the following three figures we show the separatrix of Kj and the line Ca for 
Fs = 0.0263, F^ = 0.14, flo = 0.0114, /„ = 0.2 and the initial state 4(0) = -0.3, when 
Na = 1, 3 and 35: the field values are slightly different from those used to generate 
figures IT51 and [T^ because an approximate Hamiltonian is used. 
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Figure 20 Diagram showing the curves Ca, formed by evolving the inital phase curve 
^e(O) = —0.3 through various switch times, Ta = 2iTNa/^- The dashed hues are the 
separatrixes of the resonance hamiltonian, equation 1401 at time Ta- 



For Na = 1 the line Ca is close to the initial line, /e(0) = —0.3. Slightly less than 
half of Ca lies inside the separatrix and these orbits will be transported to regions of 
larger If, and some will ionise depending upon the value of I^. 

For Na — 3 the curve Ca is more distorted; a smaller proportion of orbits lie inside 
the separatrix, but all of these are close to it and will all be transported to larger 1^ 
than the equivalent points of the previous example. 

For Na = 35, Ca has developed a complicated shape due to the motion inside the 
island. In this example a significant proportion of the orbits inside the separatrix are 
close to the horizontal line through the island centre, so will not be transported to 
regions of large /g. Clearly this structure is very sensitive to changes in Na and this 
sensitivity will be reflected in Pi. 

These figures provide a qualitative explanation for the complications seen in fig- 
ures ll8l and ll9l In particular they show why there is no simple, precise relation between 
Fa^^ and ^i"''*; they also show that the dynamics underlying the apparently simple 
resonances seen in figure^ is very complicated. 
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7 Conclusions 



In this paper we examine the behaviour of a classical hydrogen atom in parallel static 
and microwave fields, with frequencies that are low by comparison to the unperturbed 
orbital frequency. There are three main reason why the classical atom is considered. 

First, it is not necessary to make dynamical approximations in order to numerically 
integrate the classical equations of motion. The errors of the estimated ionisation prob- 
ability are determined mainly by the Monte-Carlo sampling errors and with modern 
computers these can be made acceptably small. This is in contrast to quantal cal- 
culations in which unquantifiable approximations have to made in order to solve the 
corresponding equations of motion. 

Second, within the framework of classical dynamics there is a range of easily applied 
approximations that help provide understanding of observed phenomena. The corre- 
sponding approximations are not so easy to apply to either Schrodinger's of Heisenberg's 
equations of motion. 

Finally, using techniques of analytical dynamics it is possible to construct an ap- 
proximate Hamiltonian, which provides a fairly accurate approximation to the exact 
classical dynamics, see figure O and which may be used as a basis for feasible quantal 
calculations. 

The main effects of interest here are the resonances between the microwave field 
and the Stark frequency induced by the static field. These resonances were first ob- 
served by Galvez et al (2000) and this paper also presents the first theory to describe 
these resonances qualitatively. It was shown how these resonances are responsible for 
an enhanced ionisation signal over a narrow range of static field strengths, for fixed 
microwave field amplitude and frequency. Additionally these signals disappear at par- 
ticular combinations of the two field amplitudes and recent experiments, Schultz (2003), 
Galvez et al (2004) and Schultz et al (2004), have extended the measurements of these 
"disappearance fields" . 

Since the first observations of these resonances three theoretical papers describing 
the phenomena from different perspectives have been published. Oks and Uzer (2000) 
use a Floquet analysis to derive zero-order estimates of Fs and F/i ' . Robicheaux et al 
(2002) have solved Schrodinger's equation for this problem using a split-operator method 
and have made comparisons with the classical and experimental ionisation probabilities 
for no = 39 across the j — 1 resonance with = 0.144. These calculations suggest 
that the classical and quantal values of J-s^^ are very close, but that quantal value of Pi 
is smaller than the classical value for Fs > Ts^\ In addition the time-dependence of the 
ionisation probability is described for three values of Fg and this appears to contradict 
the results summarised in figure El though it is difficult to make comparisons between 
substate averaged and unaveraged data. Ostrovsky and Horsdal-Pedersen (2003) use 
an energy shell subspace with a time-dependent electric field and weak, perpendicular 
magnetic field, with the aim of understanding oscillations seen in the experimental 
result, subsequently attributed to another cause, see Wilson et al (2004). This analysis 
is based on the same type of averaging approximation that leads to the Hamiltonian [52 
see also Born (1960, section 38), and inevitably gives zero-order estimates of Fs''^ and 

In this paper we have described the classical dynamics of this system in more detail: 
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in particular more accurate values of Fg and are determined and the properties 

of the classical resonances are described in some detail. We now list these features and 
discuss the probable consequences to the quantum mechanics. 

1) We have established that the position of the dynamical resonance, as com- 
puted by theory, is not at precisely the field J-g , at which the ionisation probability 
is largest, although Fg'''' ~ J-'s''K Further, we have shown, numerically, that the 
difference F^^^ ~ can depend on the field envelope. 

2) We have isolated the terms in the Hamiltonian that give rise to the dynamical 
resonance. This allows the computation of the dynamical resonance position using 
high-order perturbation theory; where this series converges we obtain improved esti- 
mates of the resonance position. 

This analysis is essentially the same as the corresponding quantal theory which we 
therefore expect to give the same result. Further, this suggests that the discrepancy 
noted in (QJ above will also occur in an accurate quantal calculation. 

3) Using a classical approximation, based on two stages of averaging, we have derived 
a number of conditions necessary for the dynamical resonance to affect the ionisation 
probability. These depend upon properties of the classical resonance island, the most 
significant being the island width. 

A dynamical resonance affects the ionisation probability only if it is wide enough to 
bridge the gap in phase space between the initial state and those states that ionise, 
see section 14.21 Because the ratio of these two actions is independent of the initial 
principal number, we expect a similar story in the quantal description, though quantal 
effects will inevitably blur these boundaries. This suggests that there may be cases 
where a resonance not seen in the classical ionisation probabilities will be visible in 
the quantal probabilities. 

Besides the island width, its area also plays a role in quantum mechanics; this area 
is proportional to the initial principal quantum number, ng, so we expect the reso- 
nances seen here and in current experiments to change, and possibly disappear, as no 
decreases. 

4) The resonances in the ionisation probability, not averaged over substates, are gen- 
erally very sharp, see for instance figure E| and table [3 The full-width of an iso- 
lated resonance is generally smaller than its theoretical shift produced by changing 
/„j = TO?i, as seen by comparing the data presented in tables |21 and |S1 However, 
resonances near an adiabatic boundary are significantly broadened and asymmetries 
are introduced, figure [T^ and 1161 This causes substate-averaged resonances to be far 
wider than isolated resonances, and also affects their shape. 

The classical dynamics of this process is complicated and not understood. Because the 
experimentally observed resonances behave in a similar fashion, we expect a similar 
dynamical effect in any quantal calculation, but understanding this is harder than 
understanding the classical dynamics. 

5) The shape of an isolated classical resonance and the value of Ts''\ with no sub- 
state averaging, can be affected by the field envelope, if the fields are switched on 
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sufficiently slowly, see figures El and El These changes are caused by the phase 
line representing the initial state becoming tangled as the separatrix of the resonance 
island passes through it, see figure EOl We expect this behaviour to be seen in the 
quantum dynamics provided the principal quantum number is large enough, but how 
large is not known. 

6) Classical resonances develop over a time scale that is much longer than that of 
the ionisation process operating away from resonance. This is because, at resonance, 
ionisation happens by transport around the resonance island and this is a relatively 
slow process. Further, across a resonances the ionisation time appears to reflect the 
island dynamics. For instance, at the edge of a resonance the ionisation time is longest 
because motion near the separatrix is very slow: this and other, not understood, 
features are seen in the lower panel of figure [T7I 

As with point 13 above, we should expect to see similar behaviour in a quantal calcu- 
lation, provided the principal quantum number is sufficiently large, but how large is 
not known. 
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8 Appendix: action variables 

The derivation of the required results is easiest if scaled variables are used. If a is any 
scale length and Iq an action, suitable scaled variables are 



Taking Iq to be the initial value of /„ and a the semi-major axis of the initial Kepler 
ellipse, a = In/fJ'S^ we have 



In the following these scaled units are used but for clarity the tilde is not shown. Most 
of these results are obtained using a Maple program to manipulate the series, which 
were computed to higher-order than quoted here. 

Series for Ii 

For Ii set — y so equation |S1 becomes 






(41) 
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where W = —E > 0. In the parameter range of interest fi{y) has three real roots, 
two positive Q < yi < y2 (with j/i = only when = 0) and one negative root, 

2/3 = —y-/F < 0. We may write 

h{v)^{Fy + y-){y-yi){y2-y), ymy- = 1^ {yi +y2)F = y. -2W, 

so that y- = 0(1) and, if = 0, y_ = 2W. Only the value of t/_ is needed because all 
quantities of interest can be expressed in terms of the combinations yi + 2/2 and yiy2- 
Then the series for h is 



Jyi y ]i y- \y-j 



(42) 



where 



Hk = -^ <f> dzz'^ ^h(z) and 
471- Jc 



2fc!r(3/2- fc)' 



and where h{z) = — y/{z — yi){y2 — z) has a cut on the real axis between yi and 2/2 such 
that on the real axis between these two roots on the upper branch h{x) < 0, and on the 
lower branch h > 0; for x real and x > y2, h{x) — (x — y\)[x — 1/2) and for x < yi, 
h{x) = —i\/{yi — x) (2/2 — x) . The contour C encloses the branch cut between 2/1 and 2/2 
but not the origin. 

If A: = the integrand has a pole at 2; = and contributions from the circle z = Re^^, 
as ii — » 00, so 

For A; > 1 the only contribution is from the circle at infinity, Coo 

Hk = ^£ dzz\ll-^Jl-^ = arak+i-r y\ y2'+'"^ fc>l- 

Some values are 



4:Tr Jc V z y z 



Hi = ^(2/1+^2)^-^2/12/2, i?2 = -^(2/1 +2/2) [(2/1 +2/2)^ -42/12/2] , 

^^3 = ^ [(2/1 + 2/2)^ - 42/12/2] [5(2/1 + 2/2)^ - 42/12/2] , 

^^4 = ^ (2/1 + 2/2) [(2/1 + 2/2)^ - 42/12/2] [7(2/1 + 2/2)^ - 122/12/2] • 

Perturbation expansion for 2/- 

Now we need an expression for 2/- > which is proportional to the negative root of 
/i(— 2/-/F) = 0. Because t/- = 0{1), it is expedient to define y = —zF to give the 
equation 

z^ - 2Wz'^ - 2aiFz - F^ll = 



When F = 0, z = 2W, so we put z = 2W + ZiF -\ , then perturbation theory gives 
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Thus 

+ = ^ + 4^ (^^™ - 2"?) - 1^1? (3W^^^ - 4<^?) + • • • • 

and 

y^y^ ^Ik^ ^lI^F - (Wll - 4a?) F2 + . . . . 

2VF 4VF3 161^5 V m i; ' 

Substituting these expressions for yi + ?/2 and into the above expressions for Hk 
gives 

/i = + - (2llW - 3al) - ^^"^f' - 7a?) + • • • . (43) 

2V2W 2 64VF5/2 ^ " 1^ 1024M^9/2 ^ " ^ ^ 

Series for I2 

For /2 we set i]^ = m, and equation |21 gives 

/2 = — / — y^M, /2(w) = Fu^ ~ 2Wu^ + 2a2U - 

In the parameter range of interest /2 has three real positive roots, < ui < 1*2 and 
U+/F, so we write ,f2{u) = (m+ — Fu){u — ui)(u2 ^ u) with M1M2U+ — /,?, and ui + U2 = 
(2W^-u+)/i^. Thus 



^ ^ r ^^(, _ ,^)(,^ _ /i _ !^ = ' (44) 

where and i/fe are defined above. 

The perturbation expansion for u+ is given, as before, by setting u — z/F to write 
the equation /2 = in the form — 2Wz^ + 2a2Fz — F'^I^^ = 0, giving 



This then gives 



"1 + = ^ - (W^^™ - - (3W^ - 4a^) 

and 

u^U2 = ^ + - - 4a?) F2 + . . 
and hence the expression for I2 is, 

h = 



_ i/„ _ - 3a^) - - 7al) 

2V2W 2 641^5/2 ^ " 102414^9/2 ^ " 2; 



These series for Ii and /2 now need to be inverted to give ai, a2 = 2 — ai and W 
as power series in F. The zero-order term is trivial so we substitute the series 
00 or ^ T ,00 

a,^a, + J2c,F^ a,^'-^^±^, W ^ ^ + Y^W,F\ 

k=l " " k=l 

into the series for Ii and I2 and solve for the unknown coefficients to give the energy 
and the separation constant quoted in equations 1101 and 1111 
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8.1 Angle variables 

Using the definition 9k — dS/dlk and the notation introduced in the previous section, 



Ok 



^^dy dfi/dlk ^ 1 r du df2/dlk 



By differentiating equations [S] and |51 we see that 



dy dh/dik 



= 2t:6 



Ik 



and 



du df2/dlk 



= 2TrS2k 



'vi y vA(y) 

Using the factorisation introduced in the previous section and putting 

y^Vi cos2(V'/2) + V2 sin2(V'/2), u = m cos^{x/2) + U2 sin2(x/2), 

we obtain 



(45) 



(46) 



(47) 



1 



dip {aik - Wkyi^l^) 



1 



Fy{^) 



V- 



-1/2 



dx (au + Wku{x)^ fl 



Fu{x) 



-1/2 



(48) 



where aik = dai/dlk and Wk = dW/dlk] the series expansions for both these variables 
are derived in the previous section. The first integrand can be expressed as series in 
cos mp and the second as a series in cos nx, so integration gives 



Ok = Ckiip + Pk{i>) + Ck2X + Qkix), fc = 1, 2, 



(49) 



where {Pk{x),Qk{x)) are odd 27r-periodic functions. Using the relations QHl we see that 
cii = C22 = 1 and C12 = C21 = and evaluation of these integrals, to 0(F), gives their 
Fourier series representations 



4Pi = -ai 
4P2 = -CTi 



4 - lU2I„ + 7(2/2 + Ira))F sinV' - I^afF sm2,p, 



A-ll{mn + l{2l2+ln,))F 



sin ~ I'^a\F sin 2ip, 



(50) 



4Qi = -(72 4 + /i^,(10/„ + 7(2/i+/„0)F sinx + />2'^sin2x, 



4Q2 



-(72 



4 + /3(2/„ + 7(2/i + /,„))F sinx + if^a^F sin2x. 



where ak = \/ Ik [Ik + Im)/In, k = 1, 2. When F = these expression lead to the 
formulae quoted in equation 1141 



8.2 Evaluation of dS/dF 

The generating function Si{0i,02,S.,'ri) returns to its initial value when either ^ or 77 
increases through a period, see Born (1960, page 82); here we use the notation of 
Goldstein (1980) to label generating functions. It follows that the Hamiltonian 1151 is 



44 



a periodic function of the angle variables with zero mean value. By differentiating the 
generating function S{I,^,ri) with respect to F and using the angles (V'lX); defined in 
equations 1471 we see that dS/dF may be written in the form 

i = j^f <i,>(V-2H'™ + 2„,,)(l + ^»)"'" 

1 fX / Fu\^^^'^ 
+ -— / dxiu" -2WFU-2aiF)[l , (51) 

where Wp = dW/dF and aip = dai/dF. The two integrands can be expressed in 
terms of even Fourier series in tp and x respectively where the constant terms are missing 
because of the argument given at the beginning of this section and because of the general 
result dS/dt — dSi/dt, which gives dS/dF = dSi/dF. Hence the integral |^ leads to 
the Fourier series defined in eauation llGI where the coefficients {A^, Bj.) depend upon Fg, 
Ffj, and the action variables. Because of relation 0^ it follows that no term of FdS/dF 
is independent of the angle variables. 

The Fourier series representation of dS/ dF is obtained in the same manner as that 
for Pk{ip) and Qk{x)- to 0{F), we have 

Ai = i/;iai{(/i+3/2 + 2/„) 

- i/^F (14/2 ^ 27/i/„, + 26/1/2 + 22/| + 35/2/™ + 22/^J | , 

A2 = -i/^r72{l-i^/3(4/i+6/2 + 5/™)}, As^-^FI^al 
Bi - -i/;^<72{(3/i+/2 + 2/„0 

- i/^F (22/2 ^ 35^^^^^^ ^ 26/1/2 + 2/| + 15/2/™ + 22/^) 

B2 = i/^a2{l + F/^(3/i + 7/2 + 5/„0}, B, = -j^Fl'„al 
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